Exact multilocal renormalization on the effective action : application to the random 
sine Gordon model statics and non-equilibrium dynamics 
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We extend the exact multilocal renormalization group (RG) method to study the flow of the 
effective action functional. This important physical quantity satisfies an exact RG equation which 
is then expanded in multilocal components. Integrating the nonlocal parts yields a closed exact 
RG equation for the local part, to a given order in the local part. The method is illustrated on 
the 0{N) model by straightforwardly recovering the rj exponent and scaling functions. Then it 
is applied to study the glass phase of the Cardy-Ostlund, random phase sine Gordon model near 
the glass transition temperature. The static correlations and equilibrium dynamical exponent z 
are recovered and several new results are obtained. The equilibrium two-point scaling functions 
are obtained. The nonequilibrium, finite momentum, two-time t, t' response and correlations are 
computed. They are shown to exhibit scaling forms, characterized by novel exponents Xc, as 

well as universal scaling functions that we compute. The fiuctuation dissipation ratio is found to be 
non trivial and of the form X{q''{t — t'),t/t'). Analogies and differences with pure critical models 
are discussed. 

PACS numbers: 



I. INTRODUCTION 

Recently a method was devised, the exact multilocal 
renormalization group (EMRG) [1], to obtain perturba- 
tive renormalization group equations from first princi- 
ples, in a controlled way to any order, and for arbitrary 
smooth cutoff function. It starts, as numerous previous 
exact RG studies from 
the exact Polchinski- Wilson renormalization group equa- 
tion 0, Q| for the action functional 5(0). The next 
step however consists in splitting it onto local and higher 
multilocal components |15l] . and integrating exactly all 
multilocal components in terms of the local part. This 
yields an exact and very general RG flow equation for 
the local part of the action, i.e. a function, expressed in 
an expansion in powers of the local part. 

The aim of this paper is first to develop a similar 
method using instead the effective action functional r(0). 
This is needed because T^cjj) is a very important physical 
object, both as the generating function of proper vertices, 
and related to the probability distribution of an arbitrary 
macroscopic mode (f>q . A multilocal expansion is also 
performed and yields a RG equation again in terms of 
the local part. The major advantage compared to the 
previous method Q is that one actually follows directly 
physical observables and that correlations are immedi- 
ately obtained (while in the previous method one had to 
use a second formula to compute correlations from the 
flowing action). The price to pay is a slightly more in- 
volved RG equation, but this inconvenience arises only 
at higher orders. As a simple check the 77 exponent of 
the 0{N) model will be recovered to lowest order. 

A motivation to develop such EMRG methods comes 
from disordered models. The physics of these being 
more complex than standard field theories for pure sys- 
tems, it is useful to be able to control the RG procedure. 
This is crucial for instance in the functional RG (FRG) 



which describes pinned elastic manifolds 0, 0, H^ . 
relevant for e.g . superconductors and density waves 
m in m Hi 13, and the EMRG has been apphed 
to study that problem P, . Here, and this is the sec- 
ond aim of this paper, we will study another instance of 
a glass phase, arising in the random phase sine Gordon 
model excluding vortices, as discovered by Cardy and 
Osthmd [2^. This model has been studied extensively, 
in its statics j^, 2(i, 27_, 28, _29j 3(1, _3lj and its dynamics 
m 113 , as one of the simplest but non trivial example 
of a topologically ordered glass, a continuation to two di- 
mension 21, 34] of the fixed point describing the Bragg 
glass phase in three dimension We first show that 
the present method allows to recover very simply and in a 
controlled way previous results for correlation functions 
in the statics and in the equilibrium dynamics. Next 
we obtain new results, such as the full scaling functions 
for both equilibrium and non-equilibrium dynamics. We 
obtain the corresponding exponents A and 9. We also ob- 
tain the full and non trivial behavior of the fluctuation 
dissipation ratio in the glass phase. 



The outline of the paper is as follows. First in Section 
II we derive the EMRG method for the effective action, 
and give the explicit general lowest order RG equations. 
In Section III we apply these RG equations to the pure 
0{N) model, as a test of the method. In Section IV we 
consider the Cardy Ostlund model statics. In Section V 
we study the CO model equilibrium dynamics. Section 
VI is devoted to the non-equilibrium dynamics of the 
CO model. All calculational details are contained in the 
appendices. 
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II. METHOD 
A. Exact RG method 

We want to study interacting bosonic degrees of free- 
dom described by a set of fields denoted = i/)^ where 
X is the position in space and i a general label denoting 
any quantity which will not undergo the coarse-graining 
(e.g. fields indices, spin, replica indices, additional coor- 
dinate). The problem is defined by an action functional: 



(1) 



with the initial condition Ui=o{(j)) = V{(f>), simply reflect- 
ing that the effective action equals the action when all 
fluctuations are suppressed (at ^ = where the run- 
ning propagator satisfies G/^o — 0^ from the property 
c(z,z) — 0). The above equation (PJ simply expresses 
how T{(f>) in Q depends on the final value G = Gi. The 
zero IR cutoff limit A/ = can then studied by integrat- 
ing the above equation up to / = oo. 

For actual calculations, simpler and useful choices 
read, in momentum space: 



Gl' = l(c(gV2Ag)-c(gV2A?)) 



(6) 



and by the functional integral (i.e. the partition function) 
Z = / e^"^*^"^-*. The action consists of a quadratic 
part {G^^ — G^^ is a symmetric invertible matrix) and 
V{4>) the interaction, a functional of (f). The notation : 
denotes full contractions over x,i (i.e. </> : G~^ : = 
E., <PUG-%''^y)- We wiU denote ^ / d'^x where 

d is the space dimension, and = J for integra- 
tion in Fourier. Our aim is to compute the effective action 
r(0), i.e. the generating function of proper vertices, since 
once it is known, all correlation functions are known be- 
ing simply obtained as sums of all tree diagrams drawn 
using r. For all observables to be well defined one usu- 
ally requires both an ultraviolet UV cutoff (e.g. Aq in 
momentum space) and an infrared IR cutoff (noted here 
A; = e^'Ao). For example, in a single scalar theory one 
chooses G = Gi with: 



q 



9 9 

4a2' 2A2^ 



(2) 



in Fourier. Here c(z, s) is a cutoff function which de- 
creases to zero as z ^ or s — !■ CX3 and for convenience, 
see below, we choose c{z, z) = 0. To study finite momen- 
tum observables in a massless theory one is also inter- 
ested in the zero IR cutoff limit. A; = with G = Gi=oo 
denoting c{z) — c(oo, z). 

In this paper we will use that r(0) satisfies the fol- 
lowing exact RG functional equation when the quadratic 
part G is varied (for a fixed V ((/))): 



(3) 



where the cutoff function c{x) satisfies c(0) = I and 
c(oo) = 0. With the choice c{x) = 1/(1 -I- 2x) one finds 
the massive, Pauli-Villars like, propagator: 



(7) 



with m = A^, M = Ag, where the IR cutoff mass m = mi 
is lowered from m = oo (Z = 0) to m = (^ = oo) (and 
di — > —mdm)- Whenever one needs a stronger UV cutoff 
and one may use: 



Qq 



I 



g2+r7i2''4A2^ 



(8) 



a different choice. 

The full exact RG equation Q can also be expanded 
in series of Ui as: 



which admits the graphical representation given in Fig. 1 . 

To summarize, the philosophy of the method is, in 
a sense the exact opposite of the Wilson one, since it 
amounts to start from the action with no fluctuations 
A; = Ag, and then add modes and their fluctuations un- 
til one reaches the desired theory A; <C Ag. In that limit 
one expects that the effective action reaches a fixed point 
form, given by the asymptotic solution of Q at large I. 



Derivations and more details are given in Appendix A. 
This can be used to express how the effective action 
r(0) EE Ti{(j)) of the model ^ with G = Gi depends 
on the IR cutoff A;. Indeed the following property holds: 



rii,p) = -^TrlnGi + ^^ 



4> + Ui{4>) (4) 



with Ui{4)) = Ug^{(I)) satisfies the exact flow equation: 
diUiicI,) = ^TrdiGi : (Gr' - Gr'(l + G, : ^)-') (5) 



FIG. 1; Representation of the exact RG equation @, 
The dot is the vertex Ui , the solid line a propagator Gi and 
the crossed solid line the on shell propagator dGi. The sum 
is over all one loop graphs with a factor ( — 1)''~^/2 for each p 
vertex graph represented. 
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FIG. 2: Schematic representation of the sphtting of the func- 
tional 14 vertex into local, bilocal, trilocal etc.. parts respec- 
tively (top line). Representation of the exact RG equation 
for the bilocal, trilocal, etc.. as well as local vertex (last sev- 
eral lines). Note that by definition the ERG equation for the 
bilocal part contains only exactly two feeding terms, trilocal 
three etc... The solid lines represent a propagator Gi and the 
crossed soli lines the on shell propagator dGi. Combinatorial 
factors are not represented. P here denotes the projection 
operator on the local part (denoted Pi and Pi) in the text). 



B. Multilocal expansion 

To handle the formidably complicated functional equa- 
tion Q we follow the method introduced in Q and ex- 
pand the interaction functional Ui in local, bilocal, trilo- 
cal etc.. components as: 

Uiiq})^ [ Ui{(b^)+ [ y/(0,,0„,x-2/) 



+ / Wi{(t)x,(t)y,(j)z,x,y,z) + 

J xyz 



(10) 



The local part depends only on the function (/[{(fi), 
uniquely defined from the projection operator Pi. This 
operator is fully defined in [J (see also Appendix IB|| . 



We recall here only its action on a bilocal operator 
F{cj),,cj)y,x-y), namely {PiF){cj)) = /^F(0,0,2;). It 
can be used to split an action depending only on two 
points into: 



F( 



xy 



+ / {{l-Pi)F){cj)^,d^y,x-v) 



(11) 



xy 



where, by definition, (PiF)(0, -0, z) = 5{z) F((l),tp,y), 
in such a way that the second part is properly bi-local 
(i.e. {Pi{l — Pi)F){(j)) = 0). A similar construction holds 
for higher multilocal operators . 

The idea is then to project the functional equation (jSJ 
so that the bilocal, trilocal etc.. can be expressed ex- 
actly in terms of the local part Ui only. One notices that 
there is a simplest way to do it so that the bilocal part is 
V - 0(C/2), trilocal W ^ 0{U^) etc.. This determines 
one possible splitting of the higher multilocal components 
(e.g. bilocal vs trilocal) as is represented in the Fig. |21 
and further explained in P|. This expansion is clearly 
suited to the situations where the flowing functional Ui 
becomes "small" and dominated by its local part (e.g. 
in the context of a dimensional expansion), but it has 
a more general validity, since in all cases it is an exact 
expansion in series of the local part of the full effective 
action functional. 

We now pursue the analysis exactly to order 0[Uf), 
sufficient to a number of one loop applications. Details 
are given in Appendix ^ The bilocal part is exactly 
given by: 



Vi(0i,02,a;) = ^(F,(0i,(/)2,a;) - 5{x) ^ F,((/«i, 02, Z/)) 



(12) 



with: 



Fi{c^u4>2.x) = - / di'{d\dGi.d-'){d\G'^,.d^) 
Jo 

g-ia\Gr,7«.ai-ia^Gr,7°.a^-a\Gr„.oV;,(0i);7,,(02) 



(13) 



to all orders (by definition), and the resulting exact RG 
equation for the local part of the effective action (i.e. the 
exact /3-function up to 0{U^)) is: 



l„„._n„ „ . 1 

2 



diUii^) = -aGf=O9,a,C/z(0) - - / ^Gl,^^,^kUlmGf)k„^^„AUlW (14) 



-i / d\idGf^dG^).d' fdnd\dG'f,.d'){d\Gt,.d')e-i^'-''',T''-0^~io^-Gr,r-0^-^^ 



r 



We use the following notations: di = d^, d} (resp. df) denotes derivation with respect to the first argument 
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(resp. second argument) of a function of two vectors (f>i, 
(p2, d^-G'y-d^ = J2,j G^^dld^ etc.. Also one notes in real 



space 



C'-y and Gf,i 



Gf, 



Gf = - Jl,dGf„dl". 



Note that the first line contains two one loop diagrams 
(tadpole and bubble) with one "on shell" propagators, 
and the second line represents a sum over diagrams with 
at least two loops. 



III. APPLICATION TO THE 0(N) MODEL. 

We first illustrate the method on the 0(N) model de- 
fined by ^ with 



V(0) 



92 

2 



54 

4! 



(15) 



(jjx being a N-component vector, (j)"^ = 'l2ii't'x)^ ■ The 
propagator is diagonal, and using an infrared cutoff A;, 
it reads G = Gi with: 



^1,1 J - ^ij^l 



(16) 



with Gl as in |(5J. We study this model near the di- 
mension 4, in d =^ 4 — e, and compute the effective ac- 
tion to order O(e^). For some explicit calculations, we 
will further use the form ^ with the following conve- 
nient parametrization and notation for the cut-off func- 
tion c(z): 



-t-oo 



c{z) ^ / dac(a)e-"" = / e""^ 



The condition c(0) = 1 imposes J^^ = 1. 



(17) 



A. Derivation of the /3-functions and fixed points. 

The local part of the running effective action admits 
the polynomial expansion: 



2! ^ 



94,1 

4! 



96,1 

6! 



(18) 



From power counting, it is more convenient to introduce 
the dimensionless couplings g2n.i defined from: 



92,1 = ^fg2,i 

94,1 = A;54,i 



(19) 



the coupled RG equation for g^ and g2 easily obtained by 
inserting (|18|l into H14|l as detailed in Appendix O 



Oig4,i = £34 J 5 — Ii 9 



i '94,1 



-(i[94,i 



(20) 



oig2i = 2g2i H 7 — li 94,1 5 — li g2,ig4,i 



N + 2 



^1,1' 94,1 



(21) 



with the integrals: 



f(0) _ A-2+£ 



diGl , = A^ / diGlGl (22) 

'q 



ii^J is given in ljC7|l where we show that the coefficient of 
the term proportional to ^ in (|21|l is well defined in the 
limit / ^ oo. One finds that I;^"'' — 7*^°^ is ^-independent 

and that lim;_^oo = I^^^ is universal (independent of 
c(s)) in dimension d — -i: 



m ^ L 

47r2 



dssc\s) +0{e) 



n^^^sj (2s)"^/2c'(s)(c(s)-l) 



s>0 



167r 



(23) 
2 ■ 0{e) 



where Sd is the unit sphere area divided by (27r)'' and we 
recall c'(s) < 0. Finally eq. (|20I21|) together with ^ 
lead to the fixed point values 



94 



92 



4871^ 
iV + i 



-^i'''~9l+0{e-) 



(24) 
(25) 



This fixed point describes the standard 0{N) critical sys- 
tem exactly at the critical temperature T = Tc- The ini- 
tial conditions which end up for I = oo exactly at the 
fixed point describe the critical manifold. 

Besides we obtain the correction of the critical expo- 
nent characterizing the divergence of the magnetic sus- 
ceptibility near the critical temperature from the positive 
eigenvalue A; (corresponding to the instable direction) 



diihi - 52) = ^(fe - 92) 
N + 2 

which gives correctly |35l | the exponent 7 to order e 



iV + 2 
2{N + 8) 



6 + 0(6^) 



(26) 
(27) 

(28) 



and more generally g2n,i = ff2n,iA[^ which flows to 

some fixed point values discussed below. Since 

g* ~ 0(e3) and .g*„ - 0(e") for n > 3 (see AppendixO 
for the RG equation of g^j and the free energy goj), we 
drop from now on these higher monomials and study only 



B. Computation of the 2 and 4 points proper 
vertices. 



We now compute the effective action on the critical 
manifold, up to order 0{e) for the local part, and O(e^) 
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for the bi-local part (i.e. the q dependent part), in the 
Umit of large I. Eq. 1)12(1 allows to construct the bilo- 
cal term in the effective action by inserting 1(18(1 in p3(l . 
As we restrict our analysis to order e^, we do not con- 
sider monomials higher than (0^)^ in 1(181) . and therefore 
we expand the exponential in ((13(1 to order one. Using 
the combinatorics already explained for the local part in 
Appendix O one gets 



have (for g ^ Aq): 

which coincides with the expansion of lim/^co (g) ~ 
(?^(g/Ao)~'' to order with the universal value of the r] 
exponent to this order 



^ (01 , 02, g) = ^ ^ (e'«^ - m (01 , 02 , x) 



(29) 



N ■ 



N ^ 



(3!) 



P2,X} = 



M2 



(3 



02)2) ( dl'dvGlGrgl,kf 



where we have not written terms of the form f((j)i,x) 
(i.e. which depend only on one field argument) as they 
cancel out from the effective action. To this order in e 
{0{e^)) there are no other contributions. The explicit 
expressions of Gf and diGf using 1(17(1 are given in Ap- 
pendix (|^^. This bi-local term if^ allows to treat 
the renormalization of the wave function and compute 
the exponent i] to order e^. A natural way to obtain it, 
within this method, is to compute directly the 1 particle 
irreducible (IPI) two-points function and then take the 
limit I — > oo (directly at Tc). Its local part comes from 
the quadratic contribution of ((18(1 and the bi-local part 
is the sum of G;~^(g) (Q and the quadratic contribution 

of 121 



r[S(.) 



(2; 



(30) 



0=0 

2;; 



rr(g) = Gr^(g)+ Af52,; 



N + 2 



18 



all 



In the appendix El we show that it has the form, up to 
terms of order (A//Ao)2 



rP^(.) 



Gi\q)+Af~g2,i-qHg4Mn^ 



X 



(2)ri_ 



-2 



18(47r)^ 

with the following asymptotic behaviors 



(31) 



(32) 



with a some non universal (i.e. dependent of the cutoff 
function 1(171) ) coefficient. The two-point scaling function 
X^^Hk) which is computed here (see Appendix IT^) for an 
arbitrary infrared cutoff function c{x) , is up to an addi- 
tive constant, independent of the UV cutoff I6III. F or the 
particular choice (0 one recovers the result of j3a | . 

The large argument behavior of x^"^^ (k) allows to take 
the limit Z — > 00, using the fixed point value g| ((24(1 . we 



V = V[9,] = ^^^^^ 



(34) 



in agreement with standard results |35l |. 

Let us focus on the construction of the quartic term 
in Ti((j)), obtained from the quartic contribution of l(18() 
and 1(29(1 . After combinatorial manipulations, we obtain 



r' 



quart _ ff4;Aj^ f 



4! 



K ■ 092) (093 • 094) 



_2 1 /■',,A^ + 4, 
(3!? .9. 



-9ii7^ I (( 4 )(09i ■ 092)(093 ■ 094) 

(4)/ 



093) (092 • 094))XJ (93+94) 



with X/^"* (9) defined by 



Xr'(9)- lie 



\){Gff 



(35) 



(36) 



Jqi 



and where we used the notation 

4,92,93,94 (2^)''^^'^(9i + 92 + 93 + 94). The local 
term, i.e. the first line in ((35() . contains a contri- 
bution of order which is divergent in the limit 
Z — !■ 00. Indeed, expanding it to second order gives 
gAA^I = 34,; (1 + elnA;) + 0{e^) and at first sight this 
term would lead to a divergent contribution in the limit 

I 00. However, the analysis of Xi^\q) — X^^Hl/^i) 
shows the following asymptotic behaviors 



x'^'^Hk) bk^ fc<l 



(37) 

k > 1 (38) 



with b a non universal constant. When considering the 
large I limit of the effective action, we are interested in the 
large argument behavior of x^^H^) 1(38(1 . Using the fixed 
point value gl ((24(1 , one gets that this cancels exactly the 
divergence when / — s- 00 due to the local term. Thus we 
obtain: 

i^-l™ " (nW) / ' • 094) 

4e f' , N + 4 

(( — — )(09i ■ 092)(093 ■ 094) 



(^ + 8) 4 

+ (09. ■093)(092-094))M^^^)] 



(39) 



which is independent of Aq to order 0{e^)- Note that in 
the large N limit one recovers correctly the "screened" 
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[3 7*1 four point renormalized vertex ~ eg*^ (where q is the 
transfer momentum). 

The result of this analysis is that we have constructed 
the large scale theory by obtaining directly a fixed point 
for the effective action, keeping the UV cut-off Aq finite, 
which is the relevant object for statistical physics, and 
for an arbitrary cutoff function. 



The connection between the EMRG method and the 
standard field theoretical methods in the massless scheme 
(i.e. imposing T)^'{q = 0) = 0) is more subtle here (since 
one should use / — oo strictly). 



IV. CARDY OSTLUND MODEL: STATICS 



C. Relation with field-theoretical methods 

It is interesting to make the connection with standard 
field-theoretical methods for critical phenomena. There 
one is usually interested in the limit Aq ^ oo. Note 
that in this limit H33|l diverges. It is however possible to 
define a "renormalized" effective action TE.{(f>R) which is 
well defined in that limit. 

One can first check directly on (|33|l the standard 
Callan-Symanzik (CS) "bare" RG equation 16] for the 
physical correlation function of the massless theory at the 
fixed point 

(Ao^ - r;)( hm T^M = + 0{e') (40) 

One can also connect to the CS equation for the renor- 
malized theory. One defines: 



(41) 



rfl(0) - rK\/^0) 



where Z = Z{j^,gi^i) is the so called " wave- function 
renormalization" factor such that 



+ q^ + 0{q') 



(42) 



Using 1)30(1 and noting that G; ^(q) — + 



A/ ^2^ A — <^ (0) Qj^g finds the renormal- 



0((^)^), with A = ^ 



ized mass mj^ 



J-A^C— 



|c'(0)| 



92,i) and Z 



a 



77(54,;) ln(^)). One can see that up to higher order 
terms. A/ plays the role of the renormalized mass. From 
one finds, to order {egl i,g\i) 

rriRdmnlAo In ^(^,54,;) = -di \nZ = r/lgi^i] (43) 

these derivatives being taken at fixed g^^i. This is the 
standard definition for the r](g) function. One can go 
further, define a renormalized coupling gn , e.g. through 
^^^(9 = 0) = nfpgn, with = 54 ; up to higher order 
terms, and derive the CS equations for the renormalized 
vertices. Here, we just mention one such equation 
for the "renormalized" two-point vertex function in the 
critical regime A; <C Ao but finite 



{di+ii[g^j])Tf^^\q)'^Q g/A,»l 



(44) 



obtained using the large k behavior of x^'^\k) (|32|l . We 
get again the universal value of the rj exponent form rj = 



In this section, we show how this EMRG method can 
be used to study perturbatively the Cardy Ostlund model 
[25| near its glass transition. 



A. Model, choice of propagator 

This model is a random phase Sine-Gordon model 
which can represent an XY model in a random magnetic 
field where the vortices are excluded by hand. As men- 
tioned in the introduction, the statics of this model has 
been extensively studied using various methods |0, 
mm mm The system at equihbrium is described 
by the partition function Z = / D^e^-^^^^I'^l/'^, T being 
the temperature with the hamiltonian 

H^"['l^] = l I d^x{Vdp^f- [ d^xihl cos cj), + hi sin cj),) 



(45) 

whith (f>x g] — oo,-|-cx)[ as there are no vortices, where 
= {hi, hi) is a 2d random Gaussian vector of zero 
average with fiuctuations decorrelated from site to site: 

(46) 



{hihi,)=2goAlS.,S^'\x-x') 



The quenched average over this random variable is per- 
formed by the means of replicas, which is used here as 
a simple trick to restore translational invariance and to 
organize perturbation theory. After averaging over the 
disorder, one obtains 



InZ = 

jjrcp. 



lim 



Z" - 1 



T 

ah 



a,b 

cos(C-0') 



(47) 



where a,b = 1, .., n are replica indices. We use the same 
propagator as for the 0(N) model, the Gaussian part of 
(|77|l being diagonal in replicas, one has 



T 



Gl^,^Sat-{c{q'/2Af,)^ciq'/2At)) 



(48) 



with the same decomposition of the cutoff function c{x) 
(|17|l . Notice that the hamiltonian iJ'°P possesses the 
statistical tilt symmetry (STS) j^: the last term in l|T7|) 
is invariant under the change of variable ^ -I- Ux 
which protects the diagonal (in replica space) quadratic 
term in the effective action to all orders in perturbation 
theory [Ulii- 
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B. /3- functions and fixed point 

For this model, the Fourier representation in the fields 
(|B4ll is more natural. Although only one harmonic is 
present in the starting hamiltonian H47|l . higher harmon- 
ics are generated by perturbation theory and we write 
the local interacting part of the effective action as 



El.' 



iK.i 



(49) 



where K = {Ki,...,Kn), (j) = are n- 

components vectors and one defines K.K' = KaK'^. 
The sum is over all K such that Ka are integers not all 
zero with J2a — 0. Ui{4>) is real, imposing = 5;"^' 
and the symmetry under replica indices permutation, 
which is assumed here, imposes = g^^^\ cr(K) being 
any vector obtained from K hy a permutation of the Ka- 
By inserting (|^ in (fnji (see also IjBGp in Appendix IB|) 
one obtains the RG equation for the local part to second 
order in gx- 



(1) 



47r 



P,Q,P+Q=K 



— y 



{P-Qf 



P,Q,P+Q=K 



jii J2 P Q 



(50) 



[gf = gf = gi) with P.Q = —1 After some trans- 

formations detailed in the Appendix lUl one obtains 



digi^{'^-^)9i-Bigl 
in 



(54) 



Bi = 2970 (0) / 7Ki) 



+ / (57o(i) - a7o(0))(e^=^'(^) - 1) -h 0(r)(55) 



where we used the dimensionless variable x = ccA; and 
defined: 



Gf,; = -r7^=,_p(i) 



(56) 



where the two functions d^p_(x) and 7^(2;) are given in 
(IT^ . 

As shown in the appendix, we can transform the inte- 
gral over X in (Eq. and express its cutoff dependence 
in a simple way. One finds = ^fr^iE-lJn2a) yj^j^, 

ing for T < Tc, the stable fixed point of the RG flow is 
given by 
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(57) 



with T ^ {Tc- T)/Tc and 7^ = 0.577216 the Euler con- 
stant. 



with the integrals 



Jl = A? / aCf Gf (51) 



Jl, = Ar' / (9Gr - dGr^G^GlKt 

J X 

xeJ-^ + ^''>^''T"+P-Q(^f'i (52) 

The glass transition temperature Tc below which the 
charges of minimal modulus such that i^i^-i = 
(0, .., 1, .., —1, ..,0), K1 — 2 become relevant is 



Tc = 



= 47r 



(53) 



and a small parameter t = {Tc — T)/Tc > can be 
defined, which allows to construct perturbatively the ef- 
fective action of this model H47(l in its glass phase. In- 
deed just below Tc the higher harmonics are irrelevant 
(the eigenvalues (2 — ^^-) are negative and of order 
one). Such irrelevant higher harmonics include for in- 
stance 3 replicas term Sa^&^c 
corresponding to Kl _2.i = 4- We denote gi — g]'~^ the 
coupling constant associated to Ki -i, and obtain its RG 
flow from H50|l by taking into account the 2(n— 2) possible 
fusions such that P -\- Q = Ki^-i, P, Q being themselves 
obtained by a permutation of the components of -fCi.-i 



C. Bilocal term and 2-point correlation function. 



Eq. (|B5p allows to construct the bilocal term in the ef- 
fective action to lowest order (i.e. 0(r^)) using a Fourier 
representation l|B4|l : 



Vii<l,,^,x)^Y.^i''P^e 



iK-<t>+iP-ii 



(58) 



K,P 



Just below Tcj only the charges of minimal modulus 
Kf _^ = 2 are relevant, therefore to this order the sums 
in ljB5|l are restricted to such harmonics. By inserting 
gHl into HB5|) one has 



e^^^n\K.PG^nK^,gKgP (59) 

where K,P are of the form and thus gf = gf = 

gi. Performing the integral over /' as explained in Ap- 
pendix^ we have: 



1 



KPx 



F, 



2^2 y 2^2 



7;(^)~ 



T 



(60) 
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with X = Aix. For the charges K, P we are consider- 
ing here, there are a priori 5 different cases of K.P — 
—2,— 1,0, 1,2 to consider. However we see immediately 
on the previous expression (|60|l that the charges such 
that K.P = do not contribute to the bilocal part 
of the effective action (they correspond to four rephcas 
terms .gf Ea^b^c^^d We show in Ap- 

pendixlolthat V^^'^'' takes the form, up to terms of order 
(A,/Ao)2 



V, 



^.^'^ = -A,,2(5„h.Ai+^^.P(f)) (61) 



A, 



4T4 



(62) 



where x^'^{k) behaves aymptotically at small argument 
as 



from which we extract the correlation function at the 
fixed point (up to terms of order Aq ^) 



((0,-(/)on= lim2 /(l-e*«-)[rp)]^,i 



(69) 



2Tr 



>q r \Ko) 2A2' 

2r2 In^ (|a;|Ao) + 4(1 - r + 0{t^)) In (|:e|Ao)) (70) 



which shows that the amplitude of these anomalous fluc- 
tuations in In^ (|a;|Ao) is universal (|0|)- 

We finally mention that, due to the STS, the 
connected correlation function {{4>x — (j>{Q))'^) — 
{{(j)x — (t>{0))) {{(j>x — <p{0))) is the same as in the 
pure system. 



,K.P 




K.P ^ -2 
K.P ^ -2 



fc< 1 



(63) 



V. CARDY OSTLUND: EQUILIBRIUM 
DYNAMICS 



and at large argument (relevant for the limit I —^ oo) as 



K.P 




K.P =1,2 
K.P = -1 
K.P = -2 



k > 1 (64) 



We now turn to dynamics, which, within the EMRG 
framework can be conveniently studied by introducing 
an infrared cutoff on space only, keeping the full time 
dependence. 



The large argument behavior of x^'^i^) allows to take 
the limit A; ^ of (|61() as the logarithmic divergence 
(which only exists for K ~ — P) is cancelled. We no- 
tice also that only such terms whith K = —P sur- 
vive in this limit : in particular, three replicas term as 



g«(0a-06)+»(i>b-»Ac) do not exist in the effect! 



action to order at the fixed point for A; — 0. Besides, 
by inserting the fixed point value g* (I57|l in Ai H62|l . we 
see that the cut-off dependence (encoded in the factor 
efa In 2a -J disappears in hmj^oo Ai leading to 



lim V, 



KPq 



-5K-p-z-q- In — 

lOTT Aq 



(65) 



Eq. Q together with H65|) allow to construct the bilocal 
term as 



lim r 



i- 



biloc 
I 



i4>) 



-E 

2T ^ 



(66) 



E 



I — »C30 



from which we extract the two-point 1 PI function 



r(2) .,) _ ^^r^(^) 



<p=0 
^2 



Tc{qy2Al 



(67) 



S,, + —qHn^ (68) 



47r 



An 



A. Model and propagator. 

Within this EMRG framework we want to study the 
dynamics of the model [s^ l33 |. described by a 

Langevin type equation : 



d 



V-KlUxt ^ 7 

Ot OUxt 



(71) 



where {C{x,t)) = and {C{x,t)C{x' ,t')) = 2rjTS{x - 
x')S{t — t') is the thermal noise and r/ the friction co- 
efficient. A convenient way to study the dynamics is to 
use the Martin-Siggia-Rose 's^ generating functional, on 
which perturbation theory can be done. Moreover, using 
the Ito presription, it can be readily averaged over the 
disorder. The disorder averaged generating functional 
reads 



Z[j,j] = J I?wI?rae-^["''"l+^-"+^-™ 
S[u, iu] — So[u, iii] + S'int[w, iu] 



(72) 



So[u,iu]= / iu^qtivdt + Cq )Uqt ~ VT I lUxtiUxt 
Jqt Jxt 

Sint[u, iu] = -.goAo / iuxtiuxt' cos {uxt - Uxt') 



where = dt, where in this section the initial time 
ti is sent to ti = —oo before taking A//Ao large, in order 
to describe equilibrium dynamics. 
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In our formulation the field 
ponents vector 



Uxt 

iuxt 



is now a two com- 



perturbative solution of the equation for Tiiu, iii) of the 
form ((Tn|l : 



(73) 



Ui{u, ill) 



Ui(ux,iux) 



and from 5*0 in H72|l . we compute the inverse bare prop- 
agator ^ 



iUxtFit{ux) - - 

xt ^ J xW 



(80) 

iuxt^iu'{ux)iuxt' (81) 



6{t~t'){-r]dt + cq^) 

5{t - t'){r]dt + cq^) -2r]T5{t-t') 

1 



.(g72Ag)-c(9V2An 



(74) 



where Fn(ux) and Ami^Ux) are functionals only with re- 
spect to the time dependence, i.e. functions of the "vec- 
tor" Ux = {uxt} at a given point x in space. In addition 
these will acquire an explicit time dependence, indicated 
by their t and f indices. One has the initial conditions 



By inverting this matrix we obtain the bare response and 
correlation functions 



= ^e-^^l*-'l(c(gV2Ag)-c(,V2A?)) 



(75) 



6{Uqt) 

= 9{t - t')e-?' (*-*') ^"^"2 /o A 2^ „/„2 /o A 2 



(c(gV2Ai^)-c(g72Af)) (76) 



where we have set the bare rj = 1. As we consider here 
the equilibrium dynamics of the system, the time transla- 
tion invariance (TTI) and the fluctuation dissipation the- 
orem (FDT) hold. These properties hold to all orders in 
perturbation theory and, as we will see has strong conse- 
quences on the structure of the effective action Ti{u, iu). 
This means for the dressed (i.e. exact) response and cor- 
relation functions: 



niU,^-e{t-t')-dtCiU 



(77) 

(78) 

(79) 



B. Response function and dynamical exponent. 

We will study the dynamics near the transition tem- 
perature Tc (I53|l . below which the lowest harmonic of 
the disordered pontential becomes relevant. Near Tc, we 
showed previously that the higher harmonics, although 
generated by perturbation theory, are irrelevant. 

As we considered here static disorder, the average over 
the disorder generates an effective interaction 6*1111 [u,m] 
in H72|l which is non local in time, so we expect the fric- 
tion coefficient to be renormalized by the disorder. We 
therefore construct the effective action to order one in 
r = (T — Tc)/Tc, and extract the dynamical exponent 
z from the response function. In the starting dynamical 
action H72|l the interacting part is purely local in space, so 
to order one the interacting part of the associated effec- 
tive action Ti{u, iu) will remain so. We therefore search a 



Ai^ottK"") = 25oAoCos(ut ~ Uf) 
Fi=otiu) = 



(82) 



The Fif{u) term is indeed generated by perturbation the- 
ory and is related - in the case of equilibrium dynamics 
- to A/(t/(u) by a generalized FDT relation, namely a 
Ward Identity, which can be written to lowest order: 



SUxt' 



dt'Aiu'iu) t>t' 



(83) 



where dt' acts only on the explicit time dependence (i.e. 
not on Ufi). Notice finally that terms containing higher 
powers of the field m, ie {iu)^'^'^ are of order tP+^. They 
correspond to higher cumulant of the disorder (i.e. higher 
number of replica terms in the statics). The exact RG 
equation to order one (|14|l then reads (see Appendix lEjl 



diAiu'iu) 



tit'i 



kSt'A^tAu) 



(84) 



dMu) = I k\^l,Fu{u) - / k\;i A^^,,{u) 



u(2) 



with 



(1) 
ltit[ 

(2) 
Hit', 



1 s 



^'^'M'? X 

2 duxti ' OUxt[ 



(85) 



The solution of this coupled set of equations H84|) together 
with H82() is given by 



Aitt'iu) = 2Afg,e'^''-'' cos(u:i;t - Uxt') (86) 

Uxt') t > t' 



^^''^""^ =~2Khie^^r-^'RH=?'C0s{uxt 



SUxt 

where we can check explicitly the previously mentioned 
generalized FDT relation (|83|l . Finally, as we consider 
here static disorder, the flow of gi is given by the previous 
study, the fixed point value g* being given by (|57|l . 

From Ti[u,iu\, we obtain the response function in the 
following way 



'^it-t' 



'■qtT' 



-qt' 



SiUqti^U^qt'^ 



u—iu—0 i 



(87) 



q,t,t' 
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We define 



which gives 



u=0 



(88) 



notice that in the case of equihbrium dynamics -D/tf 
Dit-t' and 'Euf = ^it-t'- One gets 



SiilqtSU-qt' 



= 5{t-t'){q^ + dt) + ^lt-t' (89) 



When considering equihbrium dynamics, the use of 
Fourier transform allows to compute this matrix element 
(|87|l in a simple way 



1 



(90) 



where 7^^J = / ^e-^'^'^nl^ and S;^ is the Fourier trans- 
form, of Y^it-t' ■ in Appendix IeI we show that it has the 
following form (up to terms of order A^/Aq) 



A2 



B, = 



9ie 



So. 2a 



2Tr 



whith the following asymptotic behaviors 



/(dyn) 



[v) ^ In S> 1 



(91) 
(92) 



(93) 
(94) 



where Odyn is a non universal constant. The large argu- 
ment behavior of x'"'^^'^^ (|94() allows to take the large I 
limit in H91|l as the logarithmic divergence is cancelled, 



lim 7^f , = — — 

i~*oc - iuj + iusB* In ■ 



B* 



lim Bi ~ e? 

I — ^oo 



(95) 
(96) 



where we have used (|57f) to compute B* which is uni- 
versal : the cut-off dependence encoded in 2a 
disappeared. On the other hand we expect that the scal- 
ing function in Fourier should read: 



za;(-^)2A 



(97) 



from scaling. If the initial model possess STS then the 
coefficient of q^ is fixed to unity. The g-independence of 
the self-energy is expected to hold only to the order in 
T that we are working at, and it should be corrected by 
higher loops. Expansion of the denominator of H95|l co- 
incides with the expansion to order t of the denominator 
of H97|l and yields the universal value of the dynamical 
exponent z 

0- 2 = 2B* = 2e''^r + 0(t2) (98) 
in agreement with previous studies. 

It is interesting in view of later applications to non- 
equilibrium dynamics, and a useful check, to compute 
this response function in the time domain. Indeed, writ- 



(2) (2)" 

ing simply the identity FJ F, 



(2) 

= I, where V\ is the 
matrix of the second functional derivatives of the effective 
action with respect ot the fields Uxt and iuxt, we obtain 
a system of closed equations for the exact response and 
correlation functions 'R-ilii and C;J(^ to order one (more 
generally, Fit{u) and A/^^/ can be bilocal in space) 



(99) 
(100) 



We remind that we have chosen the Ito presription, which 
fixes the following initial condition for the response func- 
tion 

linl7^^f_t_, = 1 

TZit^t^O (101) 

Before using these equations to study non-equilibrium 
dynamics, we show how the equation for the response 
(Eq. I99|) function allows to recover the dynamical expo- 
nent z. Using ^ together with ^ and TTI (which 



holds for equilibrium dynamics), the equation for the re- 
sponse function reads 

(9t + g2)7^^«_,, = (102) 
<y —oQ 

The limit Z — s- cx) is taken as explained in the Appendix 
IeI ljE20|) . and a way to solve this equation is simply to 
say that in the rhs, we may replace TZi^_^f,, by its bare 
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value, which is simply 9{ti — t')e~'^ (ti-t ) this term is 
already of order r. 

One expects that the response function can be written 
as: 



1™ T^lt-t' 



q'-'F^^^{r-{i-i')) (103) 



where q = q/ Aq, i = tA^, i' = t'A^, with a universal 
scaling function (up to an overall non universal scale) 
such that F'^{v) for v ^ 0. As a function it 

admits an expansion in powers of r, obtained as: 



Fp{v)^e''^{{v 



FUv)+TF'p{v)+0{r') 
l)Ei{v)e'' 



(104) 



+ 6" 



1) 



as shown in the appendix^ This is established by iden- 
tifying the direct expansion of p03|l in terms of the ar- 
gument v' = (f{t — t'): 



7^ 



F^v') + iz-2) InqiF^v') + v' F^ [v'J^Oh] 
+tF^^''{v') + 0{t^) 

with the result of solving (|l()2|l . Note that the term pro- 
portional to Inq has precisely the expected v dependence, 
a check of the calculation. Since there is an overall non- 
universal scale q Xq, F^'^{v) is defined up to a change 
in the constant p defined in the Appendix El l|E34(l . 

One can check explicitly that the scaling function in 
the time domain obtained by this second method coin- 
cides with the inverse Fourier transform of H97|l to the 
lowest order in r. The asymptotic behavior of the scal- 
ing function in the time domain is: 



{v) 



ln(l/(e'^^w)) 



OO 



v^O (106) 
(107) 



the slow time decay for z > 2, arises from the 

disorder. Notice that a similar power law tail for large 
q^t has already been obtained for the diluted Ising model 

Using the FDT we also obtain the equilibrium correla- 
tion function in the scaling regime as: 



c,v-Tr'Fc'(r(<-?')) 



c 

+ 00 



F^^'^iv) = / dwF^^^iw) 



(108) 
(109) 



We conclude this section on equilibrium dynamics by 
noticing a few interesting properties. The first one is 
an exact consequence of the scaling form (|103f) combined 
with the STS. Indeed, the STS imposes: 



lim 

t — *oo 



dt'ni- 



1 

¥ 



(110) 



Using the scaling property we showed previously, this 
symmetry (|110ll implies 



1 



dt^-'F'j.'^iqH) = -2 



10 q 

from which it follows that 



duF^^{u) = 1 (111) 



7^f,T° = / q 
( 

T 



2-KZ 



F^^^irit^t')) = 



In {t - t') 



2'Kz{t-t') 



(112) 



where we have used FDT in the last line. Note that the 
unrescaled time t appears in these formulae. Although 
the scaling form H103|) is only valid for small q we be- 
lieve that the behaviors may actually be the exact 
leading ones in the large t — t' limit, their coefficients 
being fixed (non-perturbatively) by the STS. This would 
be interesting to check numerically. 

The second property is a comparison with the so-called 
Porod's law ^J. If the form were to hold to all 
orders, the scaling functions would decay at large argu- 
ments as F'f^iv) - and F°''(w) - 1/vi. That 
yields 



w 



1 



{t - t')V'-q* 



(113) 



as in the Porod's law with d = 2 and n = 2 Hl\. Here 
this property holds to the order of our calculation 0(r). 



VI. NON-EQUILIBRIUM DYNAMICS OF THE 
CO MODEL 



Applying standard scaling arguments, we expect TZ^f, 
and C'j, to be functions of the scaling variables g^f and 
q^f where q = q/ Aq and t = Apt and z is the dynamical 
exponent. As is the case for pure systems at a critical 
point one can write from RG arguments 42] with little 
restriction: 



7^?-, = q 



n-^+^+^li] Fnir-{i^t'),i/f) (114) 



4, = rr'+"(^pj Feint ~t~'),i/i') (115) 

where the exponent 9 is defined by imposing the following 
behavior of the response scaling function Fr{v,u) when 
M ^ OO : 



Fr{v,u) ^ Fr^oc.{v) + Oiu-^) 



(116) 



This has been checked for pure systems |42l l43t l44l |45j 

and, partially for one case of a disordered system (only 
for the response function in "4^ and for the Fourier mode 
g = for both functions in 46J). It was found in all the 
pure cases that one also has 



Fc[v,u) = \-0{u ) 



(117) 
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These forms (|116I117|I yield a non trivial Fluctuation Dis- 
sipation Ratio (FDR) characterizing the violation of the 
Fluctuation Dissipation Theorem (FDT) [Uliil. It has 
been computed exactly for the spherical model in c? > 2 
|43|. using dynamical RG methods for the pure 0{N) 
model at criticality up to two loops in an e = 4 — d ex- 
pansion '45'!, and up to one loop for the critical diluted 
Ising model in a ^/e expansion 46]. 

Another standard definition for the autocorrelation ex- 
ponent Ac 0,113,1^ and for the autoresponse exponent 
Xr H is: 

Cl,^f-^Mri) (118) 
Tli~, =i^Mri) (119) 

in the limit t 00, q with t' fixed and q^t fixed, 
with (/)_R,c(0) = Cst. Assuming the behaviors (|116I117|I 
one finds the connection: 

{d-Xc)z-^ ^e-l + (2~rj)z-'^ (120) 
Xr — Ac 

which seems to hold for pure models, together with the 
inequality d/2 < Ac = A^ HHH. 

For the nonequilibrium dynamics of the CO model, we 
obtain similar scalings H114I I115|) . (77 = in this case 
because of STS) but with a different asymptotic behavior 
at large u of the scaling function Fc{v,u). As we will 
see, this has strong consequences on the FDR. Note that 
although C?- is the full correlation function, to this order 
in the r expansion it coincides with the connected one 
(which is the correct one to consider e.g. to obey FDT in 
the equilibrium regime), the difference between the two 
being of order g"^ — 0{t^). 

A. General framework 

We want to study the dynamics of the system described 
by H71I) which, at the initial time ti — 0, is in a non equi- 
librium configuration Uxti = w^, whose statistical weight 
is given by e^-f^"!""! (where Ho[uo] ^ Hco[uo])- The 
general framework to incorporate this feature in the MSR 
formalism has been developped in |42j | , and it amounts to 
describe the system in terms of the generating functional 
S[u, ill] S[u, ill] + Ho[u'^]. If the system is prepared in 
a high temperature state, with short range correlations 
— mQ^S'^{x — x'), the corresponding i?o[u°] is 

given by 

Ho[u'] = ^ l^iuD' (121) 

any addition of anharmonic terms in _ffo[u''] is irrelevant 
as long TOq 7^ 0. Moreover by power counting one has that 



mp ^ is irrelevant so that to study the leading scaling 
behavior it is sufficient to assume TOq ^ = 0, i.e. — 0. 
The effect of this nonequilibrium initial condition is then 
completely encoded in the lower bound = on the time 
integrals in the MSR functional H72() . The running bare 
response and correlation functions are given by |42l | : 

*« , = S{t - t')e-'"-''> (^e(^) - c(^)) (122) 

B. Nonequilibrium response function 

In order to compute the response function, we solve 
perturbatively the equation for TZl^, using the trick 
explained above, i.e. replacing the exact TZ'^^, in the rhs 
of H99|) by its bare value. Doing this, we obtain a per- 
turbative expansion of the exponents z (already obtained 
previously (EHl), ^ and of the scaling function Fr{v, u) in 
the same spirit as (|1U4|) . Indeed, as shown in Appendix 
IfI one has the scaling 1)1141) , in terms of the scaling vari- 
ables V — q^{i ~ t') and u = i/t' with 

Fr{v, u) = F° (t.) + tFI{v, u) (123) 

FUv,u)^F],^\v)+F],^°-^\v,u) 

e = e'^^r + 0(r2) 

which is established by comparison with the direct per- 
turbative expansion of (|114|l in powers of r: 

-R-li, = F^r{v') + (z - 2) \nq{FUv') + v'F'r{v')) 
+6ilnui^]^(u') + rF^(w',u) + 0(r2) (124) 

with v' = q^{i- i') and F]^'^{v) is given by ((Tni|l and 
Fjt,"™°'^(u, u) given in HF16P has a complicated expres- 
sion left in the Appendix (|F18|) . However its asymptotic 
behaviors, which we now focus on have remarkably sim- 
ple forms. First, in order to compare with the prediction 
for pure critical systems one is interested in the limit 
of large u, keeping v fixed. This defines Fr^oo{v) H116f) 
which, we find to be: 

Fr,oc.{v) = e-" + e^^r| - ^ Erf ^^ 
-e-'"{{l-v) In {Ave""') 

~2v{v - 1}, {| 2}, v)) } + 0(r2)(125) 

where Erf(z) is the error function and 
2F2({1, l},_j_|, 2], z^) is a generalized hypergeomet- 
ric series [SJ, IS^i IS^. This shows that the response 
function has a scaling behavior as predicted for pure 
systems at a critical point (|116|) . The small v behavior 
of Fii,oo{v) ~ 1 — e^^rlnu shows that (t>R{v) (|120|l has a 
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good limit when v ^ 0, (j)R{0) = Cst, and this gives the 
autocorrelation exponent (|119l) : 



(126) 



It is also interesting to analyze the asymptotic behavior 
in the limit of large v (and in particular v ^ u), keeping 
u fixed. This limit is relevant e.g. to study the behavior 
at fixed q, large t,t' with u — t/t' fixed. It is obtained 
from HI''16|I as explained in the Appendix|FJ The behavior 
of the response function in this limit is then given by 

lim Fr{v,u) ^ e-^' + ^Pr{u) + 0{tv~\t'') 



V — *OQ,u?ix.cd 



Pr{u) = 



u + 1 



(127) 



Notice that in the limit u ^ 1 we recover the result of 
the equilibrium dynamics H107() . This is more general as 
one can check from (|FT8)| that F'^'^''\v, u) = 0{{u-\f) 
as M — > 1. Finally, one must keep in mind that the limits 
V ^ oo and u — > oo do not commute, indeed one expects 
that a scaling function oiv/u ~ q^t' interpolates between 
these limits, left for future investigation. 

Another interesting behavior is the limit of vanishing 
momentum g = 0, the so called diffusion mode. Although 
well defined, this limit is a bit peculiar due to the pref- 
actor in the scaling function H114|) . However, the 

function Fr(v,u) behaves when u ^ in such a way 
to cancel this divergence as in H106() and leads to well 
a defined response function which has the scaling 

form 



1 



**' (t-t')(^-2)A \t' 



i) Fr (i ] (128) 



F^''{u)^Fi^^{u) + TFi''\u) 
F^iu) = 1 



i^diffi(y) ^ 26"^^ In 



1 + 



which is identified with the pcrturbative expansion of 
T?.?^'^ straighforwardly obtained from the general expres- 



sion IIFT81 



C. Nonequilibrium correlation function. 

To compute the correlation function, instead of solv- 
ing the equation for Cf^^, pU(J|l . we obtain it using the 
following formal solution for t > t' 



C?~ = lim C^^^, 

I — >-oo '■^^ 



(129) 



Jo 



where Dt-^t2 = linij^oo ^itit2 is defined in H88|l and explic- 
itly given in ljF4ll . that we expand perturbatively using 



the expression we obtained for 7^?-^, ^7^. In the Ap- 
pendix, we show that C?- has the following scaling form 
((TT3|l with 

Fciv,u)^F^{v,u) + TF^{v,u) (130) 

F^{v,u) = e-^ ^e-^^ 

and Fq{v,u) given in Appendix. Again, this is estab- 
lished by identifying the direct pcrturbative exansion of 
(ESI: 

4 = J(^5K'-) + (--2)ln(,>'^^^ 
+e\nuF^{v',u) + TF^{v',u)^ (131) 

with v' — (f'{t ~ t'), which is similar to the scaling form 
expected for pure systems at a critical point H115|) . How- 
ever, the large u behavior is different, indeed one has in 
the large u limit, keeping v fixed 



lim Fc{v,u) 



U y/U 



0{u 



-2 ru-\ 



(132) 

which decays more slowly than the predicted scaling for 
pure system at a critical point H117() . Besides, using (jF23p 
Pc,ooi^) ~ v + Oiv"^), (/>c(0) = Cst (innj), this defines the 
autocorrelation exponent Ac: 



Xc = d-- + OiT^ 



+ 0{t^) (133) 



where we have used the explicit expressions of the ex- 
ponents z (|98|l . 9 H123|) and the relation (d — Xc)z^^ — 
9 — 1/2 + {2 — ri)z^^ arising from l/^/u decay of Fq{v, u) 
()132|l . Note first that Ac is different from its trivial value 
Ac 7^ d. Besides we note that Ac 7^ Xr and finally that it 
violates the bound Ac < d/2 predicted for pure systems. 

For the correlation function it is also instructive to look 
at the asymptotic behavior v ^ 1, u fixed. As detailed 
in the Appendix IfI one has 

lim Fciv,u) F^{v,u) + -Pciu) + 0{tv-'^,t^) 



Pc{u) = Pr{u) 



(134) 



Finally we also study the correlation function in the 
limit of vanishing momentum q = 0. As mentionned 
previously for the response function this limit is a bit 
peculiar due to the q~^ prefactor in (|115|1 . The small v 
behavior of Fc{v,u) leads to the scaling form (up to a 
non universal scale) 



FS'^iu) = FS'''°{u)- 
FS^^'iu) = 1 



tF 



diffl 



in) 



F^iffi (u) = ie"^^ (4V^ + (u + 1) In (u - 1) 
-2{u~ l)ln(l + ^AI) - 21nu-^ 6 - 81n4^ 
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with the asymptotic behaviors 

F^'^iu) - 1 + re""" hi (m - 1) u 1+ (136) 
F^'^iu) - 1 + re""" ^AI u > 1 

These behaviors are such that the singularity as t' ^ 
cancels and one finds that the diffusion of the zero mode 
become anomalous at large time: 



C 



q=0 



At 



(137) 



A = 2T^ + 0{t) 



this formula being valid for t — t' ^t', the random walk 
result being recovered when z = 2. 



D. Fluctuation Dissipation Ratio. 



We now give the results for The Fluctuation Dissipa- 
tion Ratio (FDR) X^^, defined by 113 



T 



(138) 



Starting from the scaling laws that we established above, 
we can compute the FDR X^^, = X?- as a function of the 

scaling variables q^{t — t') and t/t'. As we saw previously, 
both the exponent z and the scaling function associated 
with the FDR will have an expansion in powers of r, i.e. 

-I^ = Fx{rii-i'),i) (139) 
w 

Fx{v,u)^F%{v,u) + tF},{v,u) + 0{t^) 

the expansion of z to order r being given by (|98|l . 
F^{v,u) corresponds to the Gaussian model and from 
the perturbative expansions that we obtained for 7?.?-, 

(|124|l and C?^-, H131|) one can identify (perturbatively) this 
scaling form, with 

F%{v,u) = l + e-^^ (140) 

'2^)_e^Fi^(l + e-2Tr^) 



u-l 



Fkiv,u)^ ^(1-e 

T V 

^(dF^{v,u) , u{u~l)dF^[v,u), 



dv 



du 



Inserting the formulae for and F^ obtained in the Ap- 
pendix yields the general result for Fx as a non trivial 
function of the two variables u, v. Here we only give the 
behaviors of this scaling function in the different asymp- 
totic limits studied previously. First, we note that this 
formula gives back the FDT result Fx = 1 for u = 1. 

Second, focusing on the limit u ^ 1, keeping v fixed 
one has 



Fx{v,u) = l + e~^- 



Erfi- 



(141) 



Thus in this regime X decreases below its FDT value 
XpDT — 1- Looking at this result one is tempted to 
conclude that X^^, vanishes as t/t' — > oo when q^{t — t') 
is kept fixed. In particular for q — (see below the direct 
calculation in this c ase) one finds the analogous quantity 
X^° computed in [iaiii for several models. However 
one must keep in mind that (|141|l is perturbative in r and 
the divergence of the coefficient of r could also be a sign 
of a non analyticity in r of the u = oo result. Elucidation 
of this point is left for future study. 

In the other limit that we studied previously, corre- 
sponding to u 3> 1, keeping u fixed, we obtain straight- 
forwardly the following behavior 



Fx{v,u) = 1 + e ^" 
+0{Te"v-^,T'^) 



e^^re" [u- if 
2^2 2v^ 



(142) 



This limit is relevant to study fixed q. It shows that 
there is still aging behavior in a given non zero mode, 
and appears to contradict some claims that only the 
zero mode (diffusion) exhibits interesting aging behavior. 
Note also that in this regime one has X > Xpdt, a 
feature found in other disordered models 

Finally in the limit of vanishing momentum g = 0, 
the FDR is a function of the scaling variable t/t' whose 
perturbative expansion is given by 



t 



(143) 



F^'^{u) = ™(m) + rFf ffi(M) + 0{t^) (144) 
F^'^°{u) = 2 

F^'^Hu) - 2F^'ffi(u) - 2u^^^^ - 2F^iffi(u) 



2(z-2) e 

^ ' -2- 



Zt{u — 1) T 

Using the results of previous Sections we find: 



T 



tt 



q=0 



= 2 



s'^^(\/^ + ln( 



\/m - 1 , 



cr) (145) 



where cr is a numerical constant. This constant depends 
on additive constants to respectively F^ and F^, each 
of them being nonuniversal as discussed above (see Ap- 
pendix). However, a distinct possibility is that F^^^{u) 
is universal (i.e. that the non universal parts cancel). 
Checking this can be done with the present method, 
and is left for future study. The value obtained here, 
(T = 5 — 12 In 2, may only be indicative since we did not 
keep track of all additive constants. In particular, in the 
scaling regime t 3> S> 1 , one obtains 



T 



X?r 
w 



2 + Te''"'Vu + 0{Tu°,T^) 



(146) 



Notice that taking the limit w — > (using (|F23|) ') on 
the asymptotic expression (|141|) where we have taken the 
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limit w ^ 1 before v small one recovers the same result 
(ITlHll . 

One way to understand the result 1145|) . i.e. the di- 
vergence of X^^ when t' t is to note that the same 
divergence occurs for a simple diffusion process with the 
same close times asymptotic behaviors: 

Cp'-t''/'- (147) 

nf=° ^ {t-t')^^-'^'^/' (148) 

which yields straigforwardly X^^^ A{u — 1)(2~^)/^ as 
u — > 1. Note however that to obtain the correct amplitude 
A one needs to take into account further corrections to 
Cj'jT^, specifically we note that one can rewrite H135|l as: 

Cp"^t"/'Aiu) (149) 

and that the detailed asymptotics of A{u) near u = 1 
determines the amplitude of the divergence. 

VII. CONCLUSION 

In this paper we have developped a novel EMRG 
method to perform first principle perturbative calcula- 
tions based on exact RG. Contrarily to previous works it 
is based on a multilocal expansion of the effective action 
functional. It allows to perform conveniently calculations 
with an arbitrary cutoff function in a fully controlled way 
and check explicitly the universality of the observables. 

We have tested the method on the standard 0{N) 
model. We have shown that the exponent rj to order 
O(e^) can be simply recovered within the exact RG mul- 
tilocal expansion. This is interesting since previous ap- 
proaches relied on approximations such as polynomial 
and derivative expansions, which are not needed here. 
We have also obtained several two-point scaling functions 
and explicitly checked universality. Finally, we explained 
how the method compares with more standard field theo- 
retical approaches. In a sense the present method directly 
yields the renormalized theory. 

We have applied the EMRG method to study the glass 
phase of the two dimensional random Sine-Gordon model 
(Cardy-Ostlund) near the glass transition temperature. 
We have first recovered known results for the statics 
and for the equilibrium dynamical exponent z which we 
showed to be universal. The method of derivation how- 
ever is quite different from previous ones, since it yields 
directly the self energy T,i{uj) as a scaling function of 
w/A; where A; is the infrared cutoff. We have given for 
the first time the scaling functions associated to finite 
momentum equilibrium response and correlation. 

Next we studied the out of equilibrium dynamics of 
the Cardy-Ostlund model. We obtained the two time re- 
sponse and correlations at finite momentum. These were 
found to take a scaling form and we computed analyti- 
cally the corresponding scaling functions which depend 
on two arguments v = q^{t — t') and u = t/t'. We showed 



that they exhibit aging behavior characterized by a non 
trivial fluctuation dissipation ratio X, itself a universal 
function of u, v that we obtained. We also obtained the 
off equilibrium exponents 9 and A. Interestingly we found 
that, at variance with pure systems, one must introduce 
two distinct exponents Xr and Ac for response and cor- 
relation respectively. Our study raises the question of 
whether this could be a more general property of glassy 
dynamics in disordered systems. 

Our method is promising for further RG studies of dis- 
ordered systems, as it allows to attack the problem with 
few assumptions. Other situations where it can be ap- 
plied are elastic manifolds in random media, where it can 
be used to put the so-called Functional RG on a more 
solid basis oS", '59^ . Concerning the results of the present 
paper, a numerical simulation of the Cardy Ostlund glass 
phase can be performed 60] and should provide an inter- 
esting test of the predictions of our RG calculation. In 
particular, some points require further examination, e.g. 
the asymptotic value Xoo of the FDR. This would be in- 
teresting especially in the light of the present activity on 
FDR in mean field models, and interpretations in terms 
of effective temperatures. Indeed developing real space, 
RG type methods beyond mean field remains a challenge 
in the theory of glasses. 
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APPENDIX A: EXACT RG EQUATION FOR 
THE EFFECTIVE ACTION 

Here we present a simple derivation of the exact RG 
equation satisfied by the effective action, denoted here 
rG{(f>) (and r(0) in the text), for the theory of action 
given in , when the propagator G is varied, for a fixed 
interacting functional V(0). One first introduces the gen- 
erating functional: 



ZgU) = D(t>e 



.:G- V(0)+i:0 



(Al) 



this is equivalent to the equation for Ug '■ 



S'Ug 



)tm 



(G-^ -G-\l + G 
or its equivalent form given in the text. 

Now that we have an exact equation for Fg ((/>), we can 
relate the effective action in theories with the same V((/)) 
but different G. All we need to fully determine the effec- 
tive action is an initial condition, ft is provided by the 
action itself. Indeed one has the following pcrturbative 
loop expansion: 



i.e. the partition function in presence of a set of sources 
denoted j = j^. For any variation dG of G, its variation 
dZG{j) satisfies: 



dZG{j) ^ ^\TrdG-' j 

^-\TrdG-^q4^ 
2 ojSj 



(A2) 



where dG~^ = —G~^dGG~^ and Tr denotes a trace over 
all spatial coordinates and indices x,i. Next, one in- 
troduces the generating functional Wcii) = ItlZgH) of 
connected correlations, which varies as: 



5Wg{3) 5WG{j) 



(A3) 

an exact RG equation for this quantity. From there it is 
simple to obtain the RG equation obeyed by its Legendre 
transform rG(0) = niinj(0 : j — WgU))- We will assume 
that no problem arise from the convexity condition and 
that TG{(t>) can be obtained using only the saddle point 
conditions: 



SWg 

<5Fg(0) 



{jam = <t> 



(A4) 
(A5) 



For the variation of TGi^') — (p ■ jcifp) ~ W^G(jG('/'))j this 
yields: 



5Fg(<^) = -dWGijGm 
1 1 



5] 



5j 



since the term proportional to 9g JG cancels as usual from 
the saddle point conditions HA5|I . Using HA5|I once more, 
as well as the standard relation ^^{jG{<t>)) = [^^]"\ 
gives the equation jSjl of the text. 
Writing then 



G~^ -.(fi + UGict))- ^TrlnG (A7) 



F((/)) = -irrlnG 



(A9) 



fc>i 



where F'''((/)) is the sum of all k loop IPI graphs using 
V((/)) as interaction and G as propagator. Thus if the 
initial condition for the propagator G/=o is such that all 
y^{(t>) graphs vanish when computed with Gi=o, then one 
can choose the initial condition as Ui=o{(j>) — V ((/)). This 
is the case for the choice (|6I8|) made in the text (similarly 
the initial condition for WgU) in ljA3p is the Legendre 
transform of the initial action S{(j))). 

Finally let us note that the RG equation can also be 
written as: 



d 1 



dG dG2 
2 



rrln(l + G 



S^UGi(f>) 



: (1 + G 



6^UG{^)yi 



) (AlO) 



(All) 



where the derivative in the r.h.s 



of the first equation 
is restricted to the explicit G dependence (i.e. not the 
one implicit in Ug{4>))- 



(A6) 

SWr SWr 



APPENDIX B: MULTILOCAL EXPANSION TO 

To 0(C/2) one needs only U and V ~ 0{U'^) in the 
expansion H1(J|I of U. The functional derivative reads: 



SLf 

— — 7 = S^y[dtdjU{(t)^) + I {dldjV{(l)^,(j)^,x - z) 



+dtdjV{<P,,<j)^,z - x)] + 2dldjV{<j)^,^y,x ~y) + 0{W) 

(Bl) 

using parity V{(j),'il:,—x) = V{(f>,ip,x). Inserting in 
and keeping only terms up to order 0{U'^) one finds the 
resulting RG equation: 
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+iaG?;=°(5,iaj + v^, x) + 5ia|(5Gr,yz(0, v, x) 

-6{x) I dGlVi{(^,xp,y)) + \5(:x) f 9Gf^9,9fcC/,(0)(Gf )fc™9™9,C/z(V') 

Jy ^ Jy 



where the local projection Pi operator has been ap- 
plied to obtain the first equation, and the operator 1 — Pi 
to obtain the second. This is illustrated in Fig. 2 (drop- 
ping all terms of order 0{U'^) and higher). Note that 

V{(j),4',x) = 0. The differential equation for the bilo- 
cal part Vj is linear, a general property which allows to 
solve all higher multilocal components (here V;) as a func- 
tion of the local part Ui only. The equation for Vi can 
be integrated in the forms 11 2|l . (|13|1 given in the text. 
The method is similar to 0] to which we refer for fur- 
ther details. Inserting this solution in the equation for 
Ui one obtains l|14|l in the text. We have assumed that no 
bilocal term exists in the original action. Near the fixed 
point form at large I these assumptions are not strictly 



(B2) 



(B3) 

I 

necessary, a statement which can be checked using the 
present method. 

It can be useful, in particular for the Cardy Ostlund 
model that we study in the text, to introduce a Fourier 
representation in the fields: 

V^^^^ = j d(j)dije~'^-'^-'P-'l'Vi{(l),^,x) (B4) 

Using this representation, we obtain the RG equations 
(|12I14|) in Fourier space 



Vr^" = \{Fr^ - Six) 1^ F^'-y) (B5) 

Fi^p- = - /'d/'(x.aGf,.p)(if.Gf,.p)e^^-«r'r-^+^^Gr-r ™Gr,,.p^K^p 

diU^ = ^K.dG'^-^.KU^ ~l f f {P.dG\Q){P.G^.Q)UfU^ 

^ ^ J P,Q,P+Q=K Jx 

+ i / [ P{dGf -dGf).Q f dl'iP.dGf,mP-Gf'-Q)e^''''^'^°-^''^'^ ''''^"^ (B6) 

2 JP,Q,P+Q=K Jx JO 



Where /p,Q,p+Q^^ ^ / ^^ji^SiK ^ P - Q) where N 
is the number of components of (p. In the text we have 
used V^^^^ to distinguish the Fourier series coefficients 
from the Fourier transform. 



APPENDIX C: DETAILED CALCULATIONS FOR 
THE 0{N) MODEL 



line in p4|l which reads: 



dA9o, + '-§^u' + '-fm'?\= (ci) 

\j^dGldAUi{^) - ^j^dGlGlid^d.Uim^ 



/3- function 



with implicit sums on repeated indices. Using that: 



Let us insert (|18f) into the ERG equation 114|) , keeping 
only (7o: 92 and 34 for now, and first focus on the first 



drdjUi{(j)) = g2.iSv 



9i.i 
3! 



(5,y</>2 + 2,^V^) (C2) 
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which yields: 



exactly in d = 4 



didiUi = Ng2,i + ^gl" ^ g4j0^ 

dldjUii<p,)dfdpi{<t>2) = Ngl 

-92,ig4,l((Pi+ (P2) 



(C3) 



3! 
1 



+7i^9iAiN + 4)'/'? 02 + 4(01 • hr) 



(3!)2 

Setting 01 = (j)2 in (|C3p and identifying the coefficients of 
(jy^ and (0^)^, one then easily obtains all terms in H21I2U|) 
apart from the last one, with the scaled integrals defined 
in if^ . 

Inserting now 1)18(1 into the second line of the ERG 
equation (|14|l one obtains only a correction to 32 from the 
term with the lowest number of derivatives (six). Noting 
that: 

d.djdkUi (</.) = ^ (5y + S,k^^ + 5jk<l>') (C4) 
one finds: 

iJ2dldffUMUi{<j>2) = ^-f{N + 2)01 . 02 (C5) 

i 

yielding the last term in H21|) 

-^J^dl'i^'J-gl, (C6) 



with 



//j) = / (diGf - diGf=°)dvGf,Gf,Af (C7) 



This term does not modify the fixed point value §2 to 
order e, provided it remains finite in the limit I — s- 00. A 
way to compute it is to make an integration by part to 
treat the integral over I': 



f dl'd,Gf,Grgl,Af,^ 
Jo 



(C8) 



= ^(Gr)'5L + 0(65L,5li) 

as from Eq. (|20|l 9;' (54.;' Aj,) is of order g| and where 
we have used Gf^o ^ '4'he terms we dropped are of 
order in the limit ^ ^ 00. Finally, in the large I limit 
we are left with 



dl'il'i'rgli, = ^-^Af I idiGf~diGr")iGfr + Oie^) 

(C9) 

which is already of order e^, so that the integral over x 
can be performed in d = 4 exactly. Using the decomposi- 
tion of the cutoff p7|l , we compute the following integrals 



47r2 ,1 



A-^diGf = ^ / ie--''^'/^- 



(CIO) 



Eq. l|C6p can finally be written as an integral over the 
rescaled variable x = Aix 

In the limit I ^ 00, this integral is well defined. Indeed, 
there is no UV divergence (due to the term (e~^ — 1) 
which behaves as i^) nor IR divergence due to the term 

By the same calculation one obtains the flow of the 
free energy: 

N 



dig,,= '-Af{if^~g2,i-i\'^~gli) 

^ 2 ^1,1' 94,1' 



(C12) 



= ay" / {diGf - diGr>i'Gf,Gf,Gf,M^ 



with: 



We finally obtain the flow for gg.i in H18() with the same 
kind of manipulations, and using furthermore 

9,9, ((02)3) ^ 6,5„(02)2 + 240V^02 (C13) 

9,9,9;, ((02)3) ^ 2402(^0'= + S,k^ + djkf) + 480^0^0'= 

one gets 

9,36,; = (2e - 2)ge,i - [N + U)!^''' g^jge,i 



--(3iV + 16) dl'I\^rgii, + 0{gii) (CM) 

which shows that ~ e^. Similarly there is a term 
proportional to I^^^gej in the flow equation of g^^i which 
affect the fixed point value g^ only to next order in e. 



2. Computation of the exponent 77. 

The quadratic term in (|29|) is obtained by inserting 
(|18|l in (|13(l and expanding the exponential in (|13|) to 
order one. One gets, using (|C5|) 



F,(0i,02,x)= / dl'dGt,Gf,Gf,i{Y,dldffUi{cl,^)Ui{<j>2) 



N + 2 



b2 f dl'di,Gf,Gf,Gf,rgl,Af,^ 
Jo 



(C15) 
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which is the second Hne of (|29|l . We have dropped terms 
of the form /{(j)^, x) such as 32, (54,;'/'?, h-.W^Al (resulting 
from the expansion of the exponential in l|13|l to order 0) , 
or g\i4'\, g\ (resulting from the expansion of the ex- 
ponential in (|13|l to order 1 but acting respectively with 
d^.Gfif^d^ or .Gfrf'^d^) because they do not give any 
contribution to the effective action. Indeed, the contri- 
bution of such terms to the interaction functional hii{4>) 
(Cni), will be 



Vi{(f>x,(f>v,x -y) ^ f{(j>x,x -y) - S{x -y) / fi(f)x,z) 



{/{(t^x.x -y) - 6{x -y) / ,fi(t)x,z)) 



J{cj)x,x- y) 







where we assumed parity f((j)i,x) = f{4>i, —x) (which is 
the case here) and translational invariance. To treat the 
integral over I' in IjClSp . we use an integration by part 
as in l|C8|l . one gets 



N + 2 



N + 2. 
18" 



which leads to Using (|C10|I . the last term in 

reads (forgetting for the discussion the numerical prefac- 
tor) 



H{q, Ao, A/) 
1 



x\3 



(e^^- - l)(Gr) 



(47^2)3 



(e^9^ - 1)- 



(C16) 

3 



-x^A^/2a 



-x^A^/2a 



where the integral over x is evaluated in d = 4 (as this 
term is already of order 5I For any Aq, A;, this integral 
is well defined but in the limit Ag 00, the integrand 
is not any more regularized at small x and there is a 
logarithmic divergence. We are interested in the limit 
q, Ai <^ Aq. A simple way to isolate this divergence is to 
rewrite it as 



H{q,Ao,Ai) = 



where in the second line we performed the change of vari- 
able X Aix and denoted S4 = 27r^ the unit sphere area 
in dimension d — A. Interestingly, we have (using the 
variable u — Xx^), up to terms of order X~^: 



h'{X) 



3SV 
16A 



du 



1 



-u/2a 



16A 



2a 

,-| U—OQ 



1 



-u/2a 



I _ g-'"/2a 



+ 0(A-2) = _ 



J «=0 
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where we have used c(0) — = 1, which leads to h{X) 
-^lnX + O(A-i). Finally one obtains 



H{q,Ao,Ai) 



(4^)4^ Ao 



iln^+X^^\q/A0) + O^^^ 



^-x^ I2a 



which gives (up to the factor —g\^[N + 2)/18), the last 
term in l|Sn|) . Using 1/x^ = 1/2 dtt'^e^*^^^ , one can 



compute the integral over x in x''^^ (<?) 



"9 



a.b.c J 



dt- 



-(e 4(t+°3) 



4(i + a3)' 



with ct3 — ^ + + from which we easily obtain the 
asymptotic behavior 



x(2)(^)^lng ^>1 



1 



q < 1 



as annouced in the text (|32|l . This yields a universal re- 
sult for the Tj exponent. In addition x'^^H?) gives the 
scaling function of the two point correlator x^'^\(l) = 
2q^Q{(f') where Q{y) was computed in "s^ in the par- 
ticular case of a IR "massive" cutoff function of the form 
{Tj). Although our expression is more general we have 
checked through series expansion that it coincides with 



(47r' 



/p._dhe expression given in 36] for that choice of the cutoff. 
Performing two integrations by part one can rewrite: 



\ 3 



The limit Aq ^ 00 can be taken safely in the second 
term, the UV divergence coming only from the first one 
which can be written 



X'^\q) 
1 

"3 



a,6,c>0 



a^O q q 



-gV(4a3) 



(C18) 



with C{a) = J^da'c{a'). 



-x^/2a 



-\x^/2a^3 



3. Quartic contribution to Ti((j)) 

The quartic term in H29|) is obtained by inserting p8|) 
in and expanding the exponential in to order 
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zero. One gets, using (|C3I 



Fi{(f>i,(f>2,x) 



dl'dGf,Gf,{Y,dldrfUi{^i)Ui{4,,) 



forwardly using diGf=° = !^ duc'{u) = ^. 
One has also Gf,f^ — — I). Considering specifically 
g^' ^ — gi, we first consider the possible fusion rules such 
that P + Q = : 



7V + 4 



(3!) 



(3!) 



\\2' 



1-02)') / dl'di,Gf,Gf,gi,Af,^ 
(C19) 



which is the last line in l|2H(l (here again we have dropped 
terms of the form /{(pi, x) coming from (|C3|) ). The inte- 
gral over I' in IjClOp is then treated as previously IjCSp . 
Then, when computing the Fourier transform, one ob- 
tains (jSni, with, using IClOp 



J X 



For any A;, Aq finite, this function is well defined, and we 
see that the limit Ag ^ cxd is also well defined, thus 

x!"(<l) = x'"(<l/Ai) 

the integral over x can be computed using l/x"* = 
dtte~^^ , one obtains 



dt- 



^Q'^'^ Jab Jo it + a2)^ 



.(^g 4(t + Q2) — 1) 



with a2 — l/2a + 1/26, from which we extract the fol- 
lowing asymptotic behaviors 



1 



128a2 



g< 1 



(^)(n\ 



1 
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as announced in the text (|37I38|I 



\nf g > 1 



P 
( ■ \ 
1 



-1 

V • / 



V . / 



1 

-1 

V . / 



(Dl) 



(D2) 



where . = 0, and there are 2{n — 2) different ways to 
choose P,Q like that, notice P.Q = —1. Other possi- 
ble fusions rules involve charges of higher modulus, for 
instance we could consider 



P 
I ■ \ 
1 



1 

V • / 



-1 

2 

V ■ } 



Ki.-i 



1 

V . J 



(D3) 



(D4) 



with P^ ^Q^ ^ 6. 

It is then useful to write the integrals j/"'^' and j/'^,' in 
H50|l in terms of the variables x = Aix and fi = I — I'. 
Using H56|l . and specifying to gi one has : 



APPENDIX D: DETAILED CALCULATIONS FOR 
THE CO MODEL - STATICS. 

1. /3g, -function 



^J2'9f9?iP-Q)' -in- 2)gf [ djo{iMi) (D5) 

P,Q 



The /3-function for the coupling constant g^ is ob- 
tained by inserting in (|B6|1 . This gives straight- 



and 



^Y^iP-Qf j dl'Jl.g^gf^ = {n-2)T Jid-fo{S:)-djom J df,dj^ix)Mx) - ^^{x))e 



(4-i).e^^.(^)52 (D6) 
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with X!pq' = Q p+Q=K- study the flow near 



Tf, = An, and as (jD6p is already of order g[_^, we can 
evaluate the integral over /i exactly at Tc : in particu- 
lar e^^"-^/'^-''' = 1 + 0{t). Moreover, as the integral is 
convergent, it is dominated by the vicinity of the fixed 
point n — 0. We can then substitute in l|D6ll gi-i_i by gi. 
The remaining integral over ^ is then straightforwardly 
computed by integration by parts. ljD5|l together with 
(|D6|I . integrated over and using ljD9|l . then lead in the 
limit n ^ to (|^ : 

digi = (2 - - 25f97o(0) ^ ji{x) 

/ (dloiS^) - 97o(0))(e^^^'(^) - 1) (D7) 

To compute the integrals over x in l|D7l) in the limit I 
oo at T = Tc we first quote some useful relations. Using 
the decomposition of the cut-off function H17|l , we have 



7(2a) 



in 



V(2a) 



dy 



(D8) 



a J /2a 



and the following identities: 

2x'^d^2'-ff,{x) = a7^(x) - 97^=0(2;) (D9) 

We first compute these integrals in the semi-bounded do- 
main |a;| > e and then take the limit e — > 0, in order to 
avoid problems of convergence (the limit / — > 00 does not 
introduce any problem). Let us decompose the integrals 
over X in the following way, writing as: 



570(0) / ^i{x) 



1 



(a7o(x)-57o(0)(e 



Te7oo(a;) 



1 / 97o(x)e^^^-(-) - 1 

Jx Jx 



570(0)700 (a;) 



djojx) 



970(0) (e^^''-^") - 1) 



(DIO) 



with = Using the previous formula (|D9|I for Z 

00, 2a;^9a;27oo (x) = —djfj,=o{x),x > since djoo{x) = 
0, a: > 0, together with 970 (0) = 2/Tc we are left with 
(performing the change of variable u — x"^), and denoting 

Too('^) Too('^ ) 



570(0) 



-27r 



'7^^) + ^ /'(57o(.t) - a7o(0)(e^=^-(-) - 1) 



T2 



dw(wTca„7oo(w)e^'=''~("^ -h (e^'^'^-^") - 1)) 



27r 



-I-—/ <iw(7oo(w) + u5„7oo(w)) 



-27r 

2^2 



[u(e 



Tc7oc(«) 



(Dll) 



as one recognizes total derivatives in the integrands. Us- 
ing explicitly (|D8p : 



(D12) 



7oo(^) = ^ / E^{ul{2a)) 



where Ei{z) = —Ei{—z) = J^°° e ^/z, with Ei{x) the 
exponential integral function, behaves asymptotically as 

Ei{z) ^ --fE- In z + 0{z) z<l (D13) 
Ei{z) r-. —{1 + 0(1 /zj) 2>1 (D14) 



where 7^ is the Euler constant, the limit e ^ in |ID11|) 
can be taken safely to obtain 



60 



47r 



-(7E-/„ln2a) 



(D15) 



which leads, together with (|55|l to the fixed point value 
g* (EZI). 



Bi-local term for CO model 



We compute in this section the bilocal part in the effec- 
tive action given by (|59|) . Performing in H59|) the change 
of variable V ^ ^ = I — I' and using the notations H5t)|) 
and Gf, = -T(7^(i) - 7;(x)), one gets 



KPq _ 



1 



= - I (e^^- - 1)F, 
A* 



KPx 



(D16) 



As previously, this integral is already of order gf_^ , so it 

can be evaluated at Tc, in particular e*^"'"?^'^ = 1-1- 0{e). 
Besides the integral over fi is convergent and dominated 
by = 0, so that we substitute gf_^ by gf. The remain- 
ing integral over fi is then straightforwardly computed to 
obtain 



L(g-TcK.P7l(Ai:E) _ 1) + 

(D17) 

where (ID8I1 can be written as 



pKPx ^ . 



-ii{Aix) = ^ / Ei{x^Ki/2a) - Ei{x'K'j2a) (D18) 



whith the asymptotic behaviors of Ei{z) given in (|D13I 
ID14|) . For any A;, Ag the integral over x in (|D16|I is well 
defined, but we see that in the limit Aq ^ cxd (i.e. A;, g <C 
Aq), there is a logarithmic divergence (for small x) and 
only for charges such that K.P = —2. Indeed, at small 
X, using (iDTm . -TcK.P-fiiAix) - K.P{je + lniAfx^)), 
leading to ^-'^cK.PjiiAix) ^ j^p-K.P ^ ^his imphcs that the 
limit Aq — > cxD only diverges for K.P = —2 (there is no 
problem with the large x behavior as the integrand decay 
exponentially for any couple of charges we consider here). 
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a. The case of charges K.P = 1 or 2. 



h. The case of charges K.P — —1 



For these charges, the hmit Aq ^ oo can be taken 
directly on l|D16(l . This leads to, performing the change 
of variable x — > A;x and the integral over the angular 
variable on ljD16|) 



V, 



KPq _ 2 2 



dr-(Jo(|g»-l) (D19) 







^^-K.PS^E,(r-/2a) _ ^ ^ p f E^{r^/2a)) 

J a 

where q = g/A; and Jo{z) is a Bessel function of the first 
kind. This defines the function x^'^(fc) in that case 

(e-^-^/. E^{r-/^<^) -1 + K.P [ Si(rV2a)) (D20) 



The small k behavior (the first line of lIHSl) is straighfor- 
wardly obtained as 



(D21) 



2(7b-/ ln2a) n 

aK,P = — i / u{e-^-PL _ 1 

2 J«>o 



-K.P f Ei{u/2a)) 

J a 



where we perfomed the change of variable u — r^. 

For K.P = 1 or 2, re-''-'' L Ei(rV^a) _ ^2K.p+i ^^^^^ 
r <C 1 is analytic in and using Jo(fc) ^ k^^^^ coa{k~ j) 
one finds for ^ 1 

° (D22) 
We have moreover 



dr-{M\k\r)-l) I E,{ry2a) 



2 _ 2e-"'='/2 - aP 



k^ 



(D23) 



Using (|D23|) together with (|D22|I one obtains the leading 
behavior of x^'^{k) (|D20|I in the large k limit, i.e. the 
first line of 



x''''{k)-hK.Py, 



(D24) 



bK.P = -2e2('^--/.l-2a)( / ^^-K.Pf^E.i^J _ 
Ju>0 

-2K.Pc'{0j) 

where we made the change of variable u = and used 

c'(o) = -I«. 



In that case, x^'^(^) is formally obtained as previ- 
ously (|D20|I . the small k behavior being still given by 
(ID2ip . However the large k behavior is dominated by 
the small r region and as noticed previously for r 1, 
r(e-^-^/» ^U-Vsa) _ 1) _ ^2K.p+i ^ ^-1^ ^jji^ij Yeixds to 
a logarithmic divergence in the large k limit. It can be 
obtained by computing 



drr{Jo{kr) - l){eJa Ei{rV2a) _ 



-7e+/„ In 2a 



oo 



dr 



(Mkr) - 1) 



/o r 



(D25) 



The last term in (|D20(I has the same behavior (|D23p in- 
dependently of K.P and (|D25p together with HD20|) for 
K.P = — 1 lead to the second line of 164() 



6_i = -4e''^-''"^"2a ^Q26) 



c. The case of charges K.P = —2 

As pointed out previously, there is in that case a log- 
arithmic divergence when Aq ^ 1. We isolate this diver- 
gence by writing 



V, 



KPq _ 



-1 



(qxYF, 



2 i^KPx 



(D27) 



the second term being well defined in the limit Aq —f 
oo. We focus now on the first part, using the explicit 
expression of 7/(A;a;) (|D8ll 



1 



2 fpK-Kx 



(D28) 



8T4 
-2 



e I — 1 



^ ' 2a 



El 



2aA? 



where we made the change of variable x ^ Aix. To 
anlyse the large argument behavior of 7i(A), we take the 
derivative w.r.t A 

roo 

W'(A)-^5^'y^ dxx' (D29) 



2 _ 
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where we have used E[{z) = —e ^/z. Making the change where we made the change of variable u = . The large 

k behavior is governed by the small r region in the inte- 
gral (|D33| . wherer(e2/a'Bi('-V2a)_;L_2 j^Ei{r^ /2a)) - 



of variable u = Xx^ in 7i'(A) one obtains: 
-v'iw ^ 2 2 f°° du 



using the large A behavior Ei{u/{2aX)) 
In (2a) - In {u/X) + 0{1/X) one gets 



du / -e-^e-2-f"^^(^)+0(A-2) 



-7s 



^gfgV/.(---+>"2'^)l + 0(A-2) 



(D30) 



where we have used the asymptotic behaviors l|D13ID14(l . 
This leads finally to 



{qx)'F, 



2 fpKPx 



-Sx.-pAiq^ ln( 



Ph. 
Ao' 



0(A^/Ag) 
(D31) 



which is the first term in Ht)l|l with the amplitude Ai 
given in (|62|l . 

In the second line of (|D27p . we perform the change 
of variable x — > A;x and the integral over the angular 
variable to get 



1 



{e^i-^l + -{qx)')F, 



2\ fpKPx 



1 



= -'Z'ffi ^ / dr-{M\q\r) - 1 + -q'r') 



[e'LEiir-n-) -1-2 / Ei{Ty2a)) 



(D32) 



g 27E+/^in2a^ 3^ which implies for fc ^ 1 

X^'-'^ik) ^ ± ^(Jo(|fc|r) - 1 + ifcV) + 0(1) 
-lnA: + 0(l) (D35) 
which is the last line of (|64() in the text. 



APPENDIX E: DETAILED CALCULATIONS FOR 
THE CO MODEL : EQUILIBRIUM DYNAMICS. 

1. Derivation of the RG flow. 

We restrict our analysis to order one 0{Ui), and at this 
order the RG flow reads (|14|l 



diUi(u,iu) = -dGl-fdidjUi{u,iu) 



(El) 



where Ui{u,iu) is given by H8U|) and the indices i,j for- 
mally refer to the components of the vector (j) (|73l) and 
the time dependence, i.e. di = -g^- From (|74|l . the 
matrix G^^j, has the following expression 



^-{Rf 



With these notations, we have 



2Su 



5u Su 



(E2) 



6iu 
(E3) 



where Jo{z) is a Bessel function of the first kind, from 
which we get the function x^'^(^) defined in the text for 
K.P = -2 

^^■P(k) = 4e2(7^-.Lln2a) ^ ( ( | ^| ^) _ ^ + lfc2^2) 

. n 4 



where we will often use the matrix notation for time, i.e 
u ■ V — J^utVt- Acting with this operator on Ui{u,iu), 
one gets 

^dGf=°didj{J^iuxtFit{u) - \ iu.^tiuxt' ^iw {u)) 



lUtlUf 



^ -dC[,T^.l^A,u'iu) 



(g2/^Bi(r^/2a)_^_2 / Si(rV2a)) 



(D33) 



^ dcrz^'T^Fitiu) 



the small k behavior (i.e. the second line of (|63|l in the 
text) is easily obtained 



t Ju>t[ ''out. 



X^-^(fc) - a_2A:' fc<l 



(D34) 



„2(7E-f ln2a) roc 



a-2 
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dRttTt" 



(e 



2 /^ £1(11/20) _ 



1 - 2 



duu^ 
Ei{u/2a)) 



ti>t 



(E4) 



The last term vanishes by causality since Ft{u) depends 
on with ti < t only. Identifying in (|E1|I the coefficient 
of the powers in the field iu one gets (|84ll in the text. 
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The first equation of 1)84(1 is easily solved, and it gives 
AuAu) = e'-^*^*'^ ^^'''^''-^Ai^owiu) (E5) 

where we have used C^^q — 0. From the previous study 
of the statics, one has that 



e-^'o- A^.go = Afgi + 0{gt) 



which leads together with ((E5|l to the first line in H86|l . By 
taking the functional derivative w.r.t. Uxf in the second 
line of 1)84(1 and using the same manipulations one gets 
the second line in 



2. Computation of the dynamical exponent z. 

Here we compute the self energy Yii^ 1(89(1 given by 



(E7) 



2Afgi i?r°e^«^" - d{t) / dt'R^'^e^^ 



Notice the terms proportional to S{t) in S/t (not given in 
the text ((86() for clarity) which guarantees that S/c^=o = 
0, and with the explicit expressions for the bare corre- 
lation and response functions 1(75(1 computed with the 
cut-off decompostion ((17(1 



x=0 



47r 



1 



f -\- -S^ f -I- -S^ 



(E8) 
(E9) 



After an integration by part in l(E7() and using those ex- 
plicit expressions one gets to order r 



^lui — 



— ; — lu I ate (e 

J-c Jo 



/ - 1) (ElO) 



This expression is logarithmically divergent for Aq ^ oo 
(the integrand behaves as 1/t at small t in this limit), 
and a way to isolate this divergence is to decompose this 
integral in the following way (and performing the change 
of variable t t/Af) 



V ^91. 



dte' 



'^91 ■ /alnf^^ 

-—lUJ / at[e — l)(e v '+ 2 

Tc Jq 

2gi 



1) 



Jo 



1) 



(Ell) 



where uj = w/Af and A = A^/Aq. In the first two lines 
we can take safely the limit A — *■ and we focus now on 
the divergent part of the last term 



H(A) 



T, 



c Jo 



(E12) 



(E6) Taking the derivative w.r.t. A, one has 



25/ 



I'(A) = / dt 

'c Jo 



-a 1 /.in 



H(A) 2gi iuj 

A Tc \ Jq Ja t + ^ 



dtt 



1 /a In 

-e v'^-r 



2^; iuj 



1 ^ /„ln 



— ^ / dtt / ( ^ -)e 

A ./o Ja t+^ 



(E13) 



In the integral of the second line, we can take the limit 
A — > 0, it gives 



2gi iuj 

Tc \ Jo Ja t+ ^ 



dtt 



1 Ja'n 



Tc A 7o Ja i + ^ 



0(1) 



-^!^(e/.i"(i+f)-e/oin(§)) + o(i) 
Ic A 



The last term in ((El 3(1 can be integrated by parts, to get 

Tc X Jo Jat+^ t + f 

= 5i^ + ^- /.-(^+t)+0(i) (E14) 
A Jc A 

Finally ]HI'(A) in ((El 3(1 can be written as 

H'(A)~zc.(^e/.'-2'^i + 0(l)) 

H(A) - iw(^e^"'"2° lnA + 0(l)) (E15) 
2Tc 



which gives together with the last line of l(Elip the first 
term in (jHU with the amplitude Bi = ^e-''^'"^". The 
first two lines of I(E11(I where we take the limit A ^ 
define the function x'^'^'^'^K'^) of ljnT|) : 



/(dyn) / 



r r°° 1 
(j/) = -4e-.''"'"2''| / dte*'^*(-e^"''^(*+t) - 1) 



-I-/" ^^(6*^^* - l)(ie^"'"(*+t) - 1)| 
Jo ^ 



(E16) 



The small argument behavior of 

;^(dyn)(^) is dominated 
by the large t region of the integrand (i.e. the first line of 
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(lElBll 'l. Using that (je^a _ l) ^ ^ for < > 1, 

one gets 



1 



;^(dyn)(^)__4g-/„ln2a / ^ / g^.t _ 
J a ^ 



(E17) 



which is the asymptotic behavior announced in the 
text (inS)) with the non universal amplitude cidyn — 
4e- L In I . The large v behavior of x^^^^") {v) IS gov- 
erned by the small t region of the integrand, i.e. the 
second line of 1E16II: 



idyn) 



/•I 1 

{v) - -4e- i'" ^" 2agJ„ In I / ^^(g^^.t _ _ 

-Inz^ i/>l (E18) 
which is the asymptotic behavior annouced in the text 

We show here how to take directly, in a cruder way, 
the limit / ^ oo in S;t (|E7p . Indeed, using the explicit 
expression of Cf^^ and i?fj^° (|E8|) . one has 



TJX — O 



/ (In (4A^t + 2a) - In (2a) - 21) + 0{e-'^') 



4^ t , 



- + 0(e-^') 



(E19) 



This allows to take the large I limit in at (as it is 
already of order r) 



1 



hm E^t^-^/e/o'^f^a) 

27r -/a * + 2A2 



- In (4Aot+2a) 



(E20) 



for t > 0. We then obtain directly Y^i^ in the limit I — > oo 



_*A2„f In (2a) /-oo 



lim S;^ = 

/ — >oo ZTT 



9 e 



/Jn(2a) .oo 



27r 



- ln4t+2a 



(E21) 



The small w/Aq behavior is governed by the large t region 
of the integrand, which gives 



lim S;(, 

/ — »oo 



9 e 



fjn{2a) roo ^ 



27r ./i 4t 



which gives the same result obtained by the previous 
analysis 



3. Scaling function at equilibrium 

In this section, we show how to solve the equation for 
the response function pU2|l . First, it is natural to search 



for a solution under the form TZ^ = e~'^ *Gt- Then, per- 
forming the change of variable u — t ~ ti and using the 
explicit expression ljE20p , one gets the following equation 
for Gl. 



dfG? = 45* / du 



1 



- In (4M+2a) q^u 



45* / du 



" "T 2 

1 



In (4«+2a) 



U - 



(E23) 



I 2 

2 + 



where q — q/ Aq and t ~ Apt, with the intial conditions: 



(E24) 
(E25) 



The second term in the l.h.s. is a total derivative and 
can be integrated. Performing an integration by part on 
the first term one gets 

g| = 1 + 45* (g2 f dv r due- (^"+2") e"'" 
Jo Jo 

d«e«'''e-^<'^"('^^+2'^)) (E26) 



Performing an integration by part in the integral over v 
on the first integral and peforming the change of variable 
u' = q^u in the remaining integrals one gets 

* r Jo 



f Jo 



duue^e 



+ {q t- 1)— / due 
1 Jo 



-/„l"(ff+2a) 



(E27) 
(E28) 



We now want to find the scaling function, i.e. the 
asymptotic behavior when g — > (Aq — > oo), keeping 
q^t = y fixed. In the two first lines of the above expres- 
sion, the limit g ^ can be taken safely, although the 
last term is divergent in this limit. Thus one has 

C? = 1 + AB*{{ifi - 1) / du— / due" 

* Jo 4ii 4 7o 

+ r due--f»'"(t^+'"V0(g2) (E29) 

1 Jo 

To find the asymptotic behavior of the last term we write 

(E30) 
-) (E31) 
(E32) 



1 rdue-^»'"(t^+^'^) = 
9 Jo 

1 rduie-^^'^^'P^'^^- 

r Jo 
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In the integral on the second line, we can take the limit 
g ^ by making the change of variable A = u/if', and 
the second can be done exactly. We thus have 



where we have used the expression of B* given in H95|l . 
The expression for T,itt' can be obtained in a very similar 
way : 



1 



q2 









Finally, using 
du- 



1 



4A + 2a 



1"2 



Stt' = lim Y.,tt' 



(F5) 



-AtB 



xe 



^jt-t') L - In (A^|t-t'| + f ) 

, LA2(t-t') + f 

/^ln(A^(i+i') + |)-iln(A^t+f )-iln(A^t' + f) 



1 1 



{-^E + Ei{y)~\ny) (E33) 



-6{t~t') / dh 



1 



-A§(t-ti) + f 



J„-ln(A;5|i-ti| + f) 







In (A^(t+ti)+a/2)-i In (A^t+f )- 



4u 4 
one has, up to terms of order if 

5? = 1 + - 1)-E'j(z/) + 1 - + (1 - y){\nq^ + p)) These expressions (|F4IF5|) will be very useful to deter- 



.iln(Ai^ti + f)| 



P = lE 



In - - 4 
2 



a 4A + 2a 
(E34) 



^ mine explicit expressions for 7?,^^, = lim/^oo and 
C^j, = hm;^oo by solving perturbatively I|99I1UU|) . 



which yields the scaling function given in the text. 



2. Nonequilibrium response function : detailed 
calculations. 



APPENDIX F: NONEQUILIBRIUM DYNAMICS 
OF THE CO MODEL. 

1. Some useful expressions 

To begin with, we give the explicit expression of A;tt/ 
and Ejtt' and their limiting expression when ^ — s- oo in the 
case of nonequilibrium dynamics. The general expression 
of A;tt/(u), i.e. the first line of IjESjl is still valid for non 
equilibrium dynamics. To evaluate it, we only need the 
expression of Cf^^P that we compute from H122|l using the 
same cutoff function c{z) as previously (|17|l : 



T 

^iw — ~ 



-£'°<' + '' + 2XF> + + 

Notice that the response function R^^^ has its equilib- 
rium expression. From ||E5l) and l|FlJ 

one obtains 

A/ft'(w) = e-5<^Ht -st-H't'+t-itt' A;=ott'(M) (F2) 
Using Hi^T|l one has, using T = Tc ^ in to this order: 



iim —-L^itt — T^^Wt' +^ltt' 
L — >oo L L 



-ln(A^|i-t'| + |)+ln(Ai^(t + + |) 

-iln(A^t+^)-iln(A^t' + ^)+ln^ (F3) 
and using the definition H88() . one obtains finally 



Dtv = lim But' 

I — >-oo 



(F4) 



2 

,/^-iln(Agt+f)-iln (Agt' + f) 



The starting point of our analysis is the equation l|99|l 
that we solve perturbatively by replacing, in the rhs of 
this equation, TZjf^, by its bare value. One obtains, in 
the limit I ^ oo, using (|F5|) . and in terms of the rescaled 
variables i = Apt, q = q/Ao: 



(F6) 



B* /■* 



dti-^ 



1 



(i-ii) + f 

J^ln((t+ti) + f)-iln(t+f)-iln(ti + f ) 



xe 

X ((?(ti-<')e-^'(*^-*~')-e-^' 

Let us first focus on the last term in the rhs of (|F6(I where 
we make the change of variable ti — ut and analysing the 
limit t > 1: 



-^e--(-') I duf ' ^ 



Ja-ln(l-«+^) 



L 1" In In (t+f )-i In (t) 

g/„-ln(l-«+^) 



^^al~"+^ 



2t 



xe 



J^ln(l+u+A,)-i ln(u+-i) 



(F7) 



In the integrand one can not directly take the limit t ^ oo 
because it generates a divergence of the integral when 
M — > 1. Therefore we substract the divergent term in the 
following way 



B* 



1 du 



B* 



^^al~"+^ 

^ du r 1 



/O ^ Ja ^ ^ U + 



e ' a ^ ' 2t ' 



(F8) 
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Interestingly except in the first line, one can take directly 
the limit t oo in the integrand of the two last lines 
using 



lim 



J„-ln(l-"+|f) 



2t 



J„ln(l+«+|f)-iln(«+^) 



(F9) 



The two first lines correspond to equilibrium fluctuations 
(|E23(I and their contribution to the response function has 
already been computed (|l()5|l . The last term does not 
depend any more on the cutoff function and character- 
izes the contributions coming from noncquilibrium fluc- 
tuations. The linearity of this equation suggests then to 

.^qnoncq 



look for a solution under the form TZj- — Ti-'^'p 



where Tl-^ 



i')) IHSl) and 7^?" 



w 

noncq 



e-i^(t-t')H?- with H?- determined by 



and the divergence for u ^ 1 is cured. Then all the re- 
maining integrals can be performed exactly, giving finally 



Q = B*e' 



-e-Z.^-t-fl-fOr^)) (FIO) 



We now perform exactly the same manipulations on the 
first term in the rhs of l|F6p . Performing the change of 
variable ti = ut and considering the limit i 3> 1 (keeping 
q^t,q^t' and t/t' fixed) one obtains 



B* 



1 du 
t'/t t 



't'/i t A 

B* du 



+ o{t-^) 



(Fll) 



where we have used the same trick (jF9|l as previously. 
Using (|Flip together with (|F10|) on can write HF6|I in a 
rather simple way 



g- J^ln [A(i-tt)+2a)^-if(U-i') 




(F12) 



B* 



t W 9 




(F13) 



This allows to write a close expression for the perturba- 
tive expansion of T^^noncq terms of the scaling variables 



t 



it' 



■ \VlU 



B*e-" 



— dU 



dti 



(F14) 



Unfortunately it is quite difficult to extract directly the 
asymptotic behaviors from this double integral. How- 
ever one can perform straightforward (although tedious) 
manipulations to obtain a quasi-explicit expression for 
7?.??°"'^'^. Performing the natural change of variables 

" = \/^2 — V^Ti P = + one is left with inte- 
grals over one variable 



B* 



S dt2 dti e*=-*i 



B*^~v, ^ + Q( 7,") + Q( 7,-, 

ovu iiv u — I u — 1 u 



1 rVv e^('''"i) 

"'(ITT? 



(F15) 



Pcforming further manipulations we find that one can 
write: 



e = B* + 0{t^) 



(F16) 

(F17) 
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where the logarithmic behavior determining 9 has been extracted such that the function u) has a good 

hmit for u — > oo, as will be shown below. A useful expression for this function is found as: 



^ - e- - -r^sli^^e^ (Erf/Ji^- E..,;^ 

-V^W^^e-^ (^ErfiX^-ErfijC^') + e-'" [l - v){Ei{v) -\nv - ^e) 



„ vu f 1 f V 1 \ , ^ , , 3 , V , , vu 1 , , ^ ^3 ^ vu 



-e TT — — 



+2(l-i;)e-^ln ^ + il^lnMl (FIE 

2 u — 1 J 



where Erf z is the error function, Erfiz is the imaginary 
error function: 



Erf z ^ 



Erfi z = 



dse 



dse" 



(F19) 
(F20) 



Erfz-2z/V7F 2<1 (F21) 
Erf z - 1 - e"^V(V^^) (F22) 



with Erfiz = — i Erf iz. One has the following asymptotic 
behaviors 



and 



Erfiz-2z/V7F z<l (F23) 

Erfiz - e^V(\/^^) (F24) 

and 2-F'2({l,l},{|,2},z) is a generalized hypergeometric 
series which has the following asymptotic behaviors 



2i^2({l,l},{2,2},z)^l + 0(z) 



(F25) 



,F,{{1, 1}, {-, 2}, z) + 



2 z3/2 

In(-z) 



(F26) 



Under this form, asymptotic behaviors are more easily 
obtained. Note that we have also performed numerical 
checks that (IF18|I and the starting integral l)F14(l do in- 
deed coincide. 

Note some simple formulae for the same point re- 
sponse: 

h{u)^u'^ dvFii{v,u) (F28) 
Jo 



a. Expansion at large u, v fixed. 

The asymptotic behavior of F^{v, u) is easily obtained 
in this limit. From p23|) . one has Imiu—foo Fu{v,u) — 

F ^'^ jv) + hmu^oo u ) where ^ ^^'"^(f) is given 

in (| 104(1 . On the express ions ((F16I IfT8| together with 
the asymptotic behaviors (|F21I IF23I IF25|) we see that all 
terms vanish in this limit except the following ones 

lim Fj,^°^^'^{v,u)^-Fj,^\v) 

>oo,i;rixcd 



-2v( 



'I - V^roErf ^^- e""((l - v) In (4we''^) 

,iv~^hF,i{l,l},{l'2},v))} 
which leads to p25l) in the text. 



b. Expansion at large v, u fixed. 

Although one can extract more rigorously the large v 
behavior at u fixed from the complete expression l|F16p . 
it is easier to compute it from the starting integral in 
(|F14|I . Indeed, in the large v limit, the integral will be 
dominated by the region ^2 — ~ i.e. one can replace 
in the integrand (except of course in the term e*^~*^ ) t2 
by vu/{u — 1) and ti hy v / {u ~ \): 



5*6-" f f — dt2 r dti e 



B* 1 _ 



(F29) 



2 V" 
which leads finally to 

F^™(z;,u)- 



dt2e*^ / dtie-*i 



B* 1 iV^-l)^ 
2t yfu 



-t- 0(w"3) (F30) 



We have checked that we obtain the same result by per- 
forming this expansion on (|F16|) . Finally, using the large 
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V behavior of F]^'^ H107|) and the value of B* one 
obtains 



Performing the integrals that do not involve 
one has 



Fkiv,u) ^ e-^-^^+Oiv-') (F31) 



2u2 ^ 



which gives (|127|) in the text. 



c. The limit of vanishing momentum. 

The limit g ^ is easily obtained by looking for the 
leading term in Ff>{v, u) when v = q^(t — t') — > 0. Using 
(|FT8|) together with ((1^ and (Pljl one has 



(Inu + 7b — 21n ■ 



(F32) 



This logarithmic behavior together with (|98|l cancels the 
logg divergence in (|124|l and allows to take the limit of 
vanishing momentum. We also give here the expression 

iln 
R 



of Fi"™''i(0,u), obtained from (|FTfi|l 



TTilnoncq 



(0, u) = 2e''^ In ■ 



1 



1 



(F33) 



(F37) 



r , , / dF°(v,u) ^ 
+ -(z~2)lng(^. 



2r(? 
— 

2rT 

u — 1 



2w 



-) - In 

u — 1 \u — 1 



2v 



■ IE 



dsF^ -V Fl -V, - 



u- 1 



-Fl -V Fl -V, - 



u- 1 



— e 



1 



dtl I UL2 

Jo 

1 V^T- 



this will be useful for further applications. 



3. Nonequilibrium correlation function : detailed 
calculations. 

The starting point of our analysis is the following ex- 
pression given in the text H129|l . for t > t': 



where we have used the trick (jF9p . and dropped the prime 
in v' ~ (f'{t-~t') for simplicity. A natural way to perform 
this computation is to use for F^{v,u) the decomposi- 
tion in an equilibrium and a nonequilibrium contributions 
(|123|l . Parts of (|F37|) can then be computed analytically: 



C?- = lim C'-, 

/ — !-00 '■^^ 



= 2T 



III t't\ 



(F34) 



dtl / dt^nl^Dt.t.'R-l^^ 



(F35) 



where Dt-^t2 ~ lim/^oo F)it^t2 is given in ljF4|l. that we ex- 
pand perturbatively using the expression we obtained for 
T^lt' l|f 24|l . As we did previously for the response function 
we could keep the complete cut-off dependence in HF4|I . 
However given the complexity of these manipulations and 
the experience we acquired before, we know that the only 
cutoff dependence is contained in an overall nonuniversal 
scale q Xq. For these reasons we will perform the com- 
putation using a simplified cutoff c(a) — S{a — uq) and 
we will choose uq — 2 for simplicity, i'tita can then be 
written as Hi^'4|l 



tl+t2 



ilh - t2\ + l)Vht: 



+ 0{t^) (F36) 



2Tt V 
<f w - 1 Jo 



dsFU^v\Fl^- 



u-1 



u- 1 



' 

* u — 1 



u- 1 

+ T^e""'— / dtl / dt2 











In {q^e'"')TF°{v,u) 



e / , ^ uv ^ I uv 
-4-2 Ei 



V IV 

-2 Ei 



u — 1 V u — 1 



M — 1 \M — 1 

+ ^ e -he "-!) 



(F38) 



where we have dropped the oq dependence where it turns 
out to be unimportant. 



The expressions (|F37|) together with ljF38|) allows to iden- 



30 



tify the following perturbative scaling behavior H115() 

dF°{v,u) 



One then analyses the integrals involving 

Inoncq / 



+0liiuF^iv,u) + tFI;{v,u)) + 0{t^) 



F^'"'""^{v,u) : 
2Tt V 



2v 
u-1 



(F39) 

- 7b 



uv 

4-2 -Ei 



-Ei 



u — 1 V u — 1 



M — 1 V M — 1 



u- 1 
2Tt 1 



1 - s 1 



1 - s 



Inoncq 



) 



+2 (e-"-l + t;e-'"£;i(i;)) } 



M — 1 ' S 

(F43) 
(F44) 



and the remaining integral in ljF39p where we perform 
the natural change of variable a — \/ti, /3 — 



u-1 



u — s u 



l-s 1 



u — 1 s 



I-J } 



dti 



'VhhiVh + Vh) 



f> 11 — 1 

2g2 

= 2 ^ 



da 



(^/t^~^/^)e*^+*^ 



d/3e 



a + /3 



with the exponents z and 6 given in (jHEJ and ifT^ . The The first double integral can be performed exactly 
scaling functions are universal up to a cutoff dependent 
additive constant. It was explicitly computed for the 
equilibrium response in (|E34|) . Here, we do not determine 
it and thus we can drop some multiplicative factors of 
momentum in the log q term. 



T, e 



da 



(^Erfi 



2q2 



u-1 



-Erfi 



u- 1 



a. Expansion at large u, v fixed. 

First, one has 

= e-''-e-"^ - ^^ + 0(u-2) (F40) 

We now focus on the asymptotic behavior of Fq{v,u) 
for large u, keeping v fixed. Using the small argument 
behavior of Ei(z) In z + 7^; + 0{z) one has for the first 
line of Hl''39|l in this limit 



X Erfi 



M- 1 



'J-ErfiV^+0(u-i) (F45) 
V u 



And we expand the second one in the following way 



e 



da 



2e^''t 



i^^(^)-lnf^ 
u — 1 V u — 1 



■ IE 



E 



O(u-i) 



q- 

V 

u 



dfie" 



dae-U^^- 
\2au 



a + P 



n/2 



1^(1 + 0(a)) 



(F41) 



Again using the small argument behavior of Ei{z), one 
has 



n>2 

2e''^TcT _„v da 
' 7^ e 



-Ei 



u — 1 V u — 1 



+2 {er" - 1 + ue 
„ In u 



uv 

-4-2 Ei , 

?i — 1 \M — 1 



e""- Inu + 0(u 



(F46) 



Finally, (|F40IF42IF45IF46p lead to the asymptotic fol- 
lowing form for Fc{v, u) in the limit u 00, v fixed, 

lim i-cl^^,"") = hr — ^ 0[u ,TU ,t) 



+ 0(u-i) 



(F42) 



i^c.^(u) = e^^e"''\/^Erfi/^ 



(F47) 
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notice that the subdominant terms in hiu/u cancel be- 
tween (|l''42ll and (|l''46|l so that the leading corrections are 
of order u^^. (|F47|I gives the asymptotic behavior given 
in the text 11321). 



b. Expansion at large v, u fixed. 



c. The limit of vanishing momentum. 

To obtain the limit of vanishing momentum q ^ 
of the correlation function, we look at the behavior of 
Fc{v,u) when u — > 0, up to order 0(v) terms (due to 
the q^^ prefactor in (jF39|) '). This is done in the following 

way 



In this limit, the terms in the four first lines of ljF39|) 
decay exponentially in this limit. The fifth line however 
(which corresponds to the equilibrium contribution) de- 
cays like a power law. Indeed, using the large v behavior 
of Ei{v) - e'"{l/v + + 0{v-^)) one has 



+ 0(w-2) (F48) 



We now analyse the behavior of the terms involv- 
^ in ^EM- Usir 

(w, u) HF30() . one has 



ing pjl^°^°'i in ||1''39|| . Using the large v behavior of 

Filnoncq 
; 



u-1 



ds 





1 - s 

u-1 
1 



u- 1 



V F, 



Inoncq 



R 



1-S 1 

u ~ 1 s 



Inoncq 



R 



u — 1 s J ) 



(F49) 



u-l 



ds 







(^,) {u-lf (y/TT^-l)^ 

2 /. /77TZ _ 1^2 . 



Notice first on this expression that we are left with con- 
vergent integrals over s. Moreover, in the large v limit, 
due to the exponential prefactors the first term decays 
also exponentially (for u > 1), and the second one is 
dominated by s = 1, which leads to a power law decay 

v{u-l) Jo 



We are now left with the double integral in ljF39p . which 
is dominated - also due to the exponential prefactor - by 
ti ~ vu/ {u — 1) and t2 ^ v/{u—l). Therefore to get the 
leading behavior, we substitute ti and t2 by these val- 
ues in the integrand (except of course in the exponential 
e*i+*^). This yields 



2 



dti 







dU 



t2e 



tl+t2 







2 Ju V 



dtie*- 



-0( 



1 



L Ju V f 



(F50) 



which together with the other term in w ^ ljF48p yields 
(|134|l in the text. 



26-^-6-"^ ( Ei{^) - In f ^ 
' u~ 1 \ u — 1 



IE (F51) 



!'(" + !) / , ^ uv ^ uv 

-4-2 Ei 



V ^ V 

-2 Ei 



u — 1 \u — 1 



u — I \u — I 
+ 2(e~^ + e"^) 



u-1 

,2 



(6-21nwe''^ - u Init + (it + 1) In (u - 1) 



Then using the expression of F^"°"°"^(0, m) (|F33|1 . one 



has 



1 - s 1 



u — 1 s 



,0/1 ^ \ TTilnoncq ^ ^ ^ 
' V \ Fr, 



u-1 



2v f 1 



u 

+0{v^ 
= 4e'^^ - 



u — 
1 
s 



hi)} 



(F52) 



nonoq , I ^^^noncq (^Q,^) } 



Snii^ + mi- 



u-1 Jo 

- -(V^ - (u - 1) In (1 + ^) - 2 In 4) + 0{v^) 

u — 1 ^u 

To treat the double integral in the last line of ljF39p we 
come back to the variables t,t',q 



dti / ai2 , , , — 1=7 

irre^-e-« ) / dt, / dt2 ' 



Jo 



2 



(F53) 

Under this form, the limit g — > is very simply obtained: 



lim — -re "^e "-i 
q^o 2 g2 



dt. 



1. 



Tre'"' / dti / dt 



dto 

VhhiVh + Vh) 



2TTe"'i' {{u - 1) In (1 -f V") -^^^^ 



(F54) 
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Finally l|F5Tl lF52l IF54|| . together with ||F39)| and the 

complete expression of the correlation function C?~ (|131|) 
lead to 

= 2ri' (l + r - i^(l„ (t - (■■) + Te) + 9 In i 

^^diffi (y) ^ 1 (^4^ + (u + 1) In (u - 1) 
-2(m - 1) In (1 + V^) - 2 In u + 6 - 8 In 4) (F55) 
where we have used Fq{v,u) — v/{u — I) + O(w^), 



d^F^{v,u) = l/(it - 1) + 0(v) and v/{u - 1) = g^f' ; 
this gives the scaling form p35|) given in the text. 
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